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INVESTIGATIONS  TO  DEFINE  ACCEPTABILITY  TOLERANCE 
RANGES  IN  VARIOUS  REGIONS  OF  COLOR  SPACE 

I.  INTRODUCTION 


The  establishment  of  acceptability  tolerance  limits  for  purchased  goods 
is  a  problem  as  old  as  the  practice  of  colorimetry.  Basing  the  tolerance  on 
a  perceptibility  criterion  has  become  discredited,  as  it  is  now  generally 
realized  that  perceptibility  and  acceptability  are  not  the  same.  A  method 
of  setting  up  tolerance  limits  that  really  represent  the  !•  .ormed  judgment 
of  experienced  inspectors  is  greatly  to  be  desired.  This  need  becomes  all 
the  more  cogent  as  instrumental  color  measurement  gradually  displaces  visual 
methods. 

The  United  States  Army  has  been  accepting  shipments  of  textiles  to  be 
made  into  army  uniforms  on  a  visual  basis.  The  inspector  is  provided  with  a 
standard  sample  and  eight  limit  samples  representing  maximum  allowable  de¬ 
viations  from  the  standard  in  each  of  eight  directions  in  color  space.  Armed 
with  these  nine  samples,  the  inspector  examines  the  submitted  sample  and 
issues  a  pass  or  fail  verdict.  This  visual  method  has  proved  troublesome, 
as  evidenced  by  frequent  disagreements  between  the  Army  and  the  textile  mills 
submitting  the  samples.  The  Army  therefore  wishes  to  change  its  acceptance 
procedures  to  one  using  instrumental  measurement  and  valid  tolerance  limits 
for  acceptance. 

The  object  of  this  investigation  was  to  develop  a  procedure  for  setting 
up  these  valid  acceptance  limits  for  use  with  instrumental  measurement.  The 
Army  submitted  the  standards  and  visually  selected  limit  samples  for  twelve 
colors  of  long-term  interest.  In  addition,  some  200-300  samples  were  sub¬ 
mitted  for  each  of  three  other  shades.  These  samples  were  routine  submis¬ 
sions  for  acceptance  or  rejection  received  from  various  textile  mills.  We 
were  not  requested  to  set  up  tolerances  for  any  of  these  standards;  we  were 
rather  asked  to  develop  a  modus  operand!  that  the  Army  could  use  to  set  up 
such  tolerances  for  any  sample  in  the  future. 

II.  TESTS  ON  STANDARD  AND  LIMIT  SAMPLES 

A.  Relationships  between  standard  and  limit  samples  in  color  space. 

Ideally,  the  standard  and  the  eight  limit  samples  for  any  color  should 
form  a  three-dimensional  figure  in  color  space,  with  the  standard  somewhere  in  the 
interior  of  the  figure  and  all  the  limit  samples  on  the  boundary.  We  wished  to 
see  if  this  was  indeed  the  case  for  limit  standards  that  had  been  chosen  visually 
and  been  used  successfully  for  several  years  in  Government  procurement. 

We  obtained  measurements  of  all  the  standard  and  limit  samples  on  a  Diano- 
Hardy  spectrophotometer,  together  with  calculated  CIELAB  coordinates  L*,  a*  and 
b*.  We  then  wrote  a  computer  program  that  created  a  three-dimensional  figure  In 
CIELAB  color  space  for  each  of  the  twelve  colors.  This  figure  was  constructed 
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from  the  standard  and  the  eight  limit  samples  all  taken  as  a  group.  The  pro¬ 
gram  in  effect  created  all  possible  tetrahedra  that  could  be  constructed  from 
all  the  nine  samples  taken  four  at  a  time,  a. id  then  put  these  tetrahedra  to¬ 
gether  to  form  the  three-dimensional  figure  The  program  then  determined 
which  of  the  nine  samples  (the  standard  and  the  eight  limit  samples)  lay 
inside  the  figure  and  which  lay  on  the  boundary.  Appendix  A  describes  the 
program  and  presents  the  mathematical  details,  as  well  as  lists  of  the  computer 
program  itself  (program  RAMSLEY), 

Table  1  presents  the  computer  printout  from  this  program.  For  each  of 
the  twelve  shades  the  printout  shows  the  designations  of  the  limit  samples, 
the  CIELAB  coordinates  of  the  limit  samples  as  well  as  the  standard,  and  an 
indication  of  whether  or  not  each  of  the  standards  and  limit  samples  was  on 
the  boundary  of  the  three-dimensional  figure.  We  see  that  in  only  three  of 
the  twelve  cases  did  the  standard  lie  inside  the  figure;  in  the  other  nine 
cases  it  was  cn  the  boundary.  Furthermore,  in  ten  out  of  the  twelve  cases 
one  or  more  of  the  limit  samples  fell  inside  the  figure;  in  only  two  sets 
were  all  the  limit  samples  on  the  boundary. 

Figures  1  and  2  show  CIELAB  color  space  plots  of  the  standard  and  the 
limit,  samples  for  Olive  Drab  7,  untreated  cotton  duck,  8.25  oz.  Figure  1 
shows  a  plot  of  b*  against  a*.  The  isohue  line  indicated  on  the  figure 
is  the  line  joining  the  standard  with  the  origin;  the  isochroma  line  is  a 
segment  of  a  circle  drawn  through  the  standard  around  the  origin.  The 
letters  shown  are  abbreviations  for  the  designations  of  the  limit  samples; 
their  meaning  can  be  inferred  from  the  designations  listed  for  this  sample 
in  Table  1.  Although  they  are  useful  to  visual  graders,  we  see  little 
correlation  between  the  color  names  applied  to  some  of  the  limit  samples 
and  their  position  on  this  diagram.  In  Figure  2,  which  is  a  plot  of  L* 
against  a*,  we  see  tnat  there  is  good  correlation  between  the  thin  or 
full  designation  and  the  lightness  of  the  sample;  the  thin  samples  are 
uniformly  lighter  than  standard  and  the  full  samples  darker.  A  glance 
at  the  L*  values  in  Table  1,  however,  shows  that  this  correlation  does 
not  strictly  hold  in  all  cases;  there  are  cases  where  the  thin  samples 
arc  darker  than  standard  or  the  full  samples  lighter. 

Summarizing  the  results  of  this  section,  we  see  that  the  standard 
and  the  limit  samples  do  not,  in  many  cases,  have  the  expected  orientation 
of  standard  in  the  middle  and  limit  samples  on  the  outside.  Also,  the  de¬ 
scriptive  designations  of  the  limit  samples  do  not  always  show  a  correla¬ 
tion  with  their  plotted  points  in  CIELAB  color  space.  These  factors 
probably  reflect  the  experience  of  dyers,  their  practices  in  matching 
shades,  and  their  terminology. 

B.  Determination  of  minimum  ellipsoid  containing  standard  and 
limit  samples. 

1.  Hon-tilcing  ellipsoid,  unconstrained  axis  angle.  In  order  to 
learn  more  about  the  orientation  of  the  standard  and  limit  samples  in  CIELAB 
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Tabl.  1.  Designal  ions,  CIEL4B  Coordinates,  and  Geometric 
Orientations  o£  Standards  and  Limit  Samples  (D75,  CIE  1SJ1) 


OLIVE  DRAB  7t  UNTREATED  COTTON  DUCK,  8.25  OZ 


SAMPLE 

L 

A 

B 

POSITION 

STANDARD 

THIN  STANDARD 
THIN  YELLOW 

THIN  GREEN 

THIN  RED 

FULL  STANDARD 
FULL  YELLOW 

FULL  GREEN 

FULL  RED 

34.12 

34.45 

36.05 

36.25 

36.56 

32.65 

33.19 

33.53 

32.63 

-2.71 

-2.64 

-3.15 

-2.70 

-2.22 

-2.34 

-2.49 

-2.73 

-1.70 

10.85 

11.07 

11.75 
11.04 

10.75 
10.65 
10.04 
10.82 
10.13 

ON  BOUNDARY 

NbT  ON  BOUNDARY 
ON  BOUNDARY 

ON  BOUNDARY 

ON  BOUNDARY 

ON  BOUNDARY 

ON  BOUNDARY 

ON  BOUNDARY 

ON  BOUNDARY 

tan  46»  COTTON  POPLIN,  4  OZ. 


SAMPLE 


L  A  B 


STANDARD 
THIN  STANDARD 
THIN  RED 
THIN  BLUE 
FULL  STANDARD 
FULL  YELLOW 
FULL  RED 
FULL  nLUE 


68.90 

1.03 

71.13 

.69 

67.88 

1.15 

69.57 

1.20 

67.72 

.76 

68.40 

.31 

67.10 

1 . 30 

69.23 

.81 

13.69 

13.47 

13.21 

12.13 

13.20 

13.03 

12.86 

13.09 


POSITION 

ON  BOUNDARY 
ON  BOUNDARY 
NOT  ON  BOUNDARY 
ON  BOUNDARY 
ON  BOUNDARY 
ON  BOUNDARY 
ON  BOUNDARY 
NOT  ON  BOUNDARY 


TAN  M-l t  CL.  POLY/WOOL  TROP.  9  OZ. 


SAMPLE 


L  A  B 


STANDARD 
THIN  STANDARD 
THIN  YELLOW 
thin  GREEN 
thin  RED 
FULL  STANDARD 
FULL  YELLOW 
FULL  GREEN 
FULL  RED 


54.86 

55.28 

55.99 

54.67 

54.74 

54.76 

54.42 

54.65 

54.19 


2.96 

15.73 

2.77 

15.75 

2.62 

15.82 

2.75 

15.30 

3.15 

15.62 

2.85 

15.99 

2.92 

16.06 

2.41 

15.37 

3.17 

15.  BS 

POSITION 

MOT  ON  BOUNDARY 
NOT  ON  BOUNDARY 
ON  BOUNDARY 
ON  BOUNDARY 
ON  BOUNDARY 
NOT  ON  BOUNDARY 
ON  BOUNDARY 
ON  BOUNDARY 
ON  BOUNDARY 


ARMY  GREEN 

SAMPLE 

STANDARD 
THIN  STANDARD 
THIN  BLUE 
THIN  GREEN 
THIN  RED 
FULL  STANDARD 
FULL  BLUE 
FULL  GREEN 
FULL  RED 


r  WL  GAB .  r  130Z. 

L  A 

26.49  -3.3o 

26.60  -3.39 

26.51  -3.53 

26.20  -3.61 

26.53  -3.11 

25.87  -3.15 

25.83  -3.38 

26.13  -3. SI 

26.34  -3.09 


B  POSITION 

1.34  ON  BOUNDARY 

1.12  ON  BOUNDARY 

1.05  ON  BOUNDARY 

1,07  ON  BOUNDARY 

1.52  ON  BOUNDARY 

1.30  ON  BOUNDARY 

.88  ON  BOUNDARY 

1,42  ON  BOUNDARY 

1.6?  ON  BOUNDARY 
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Tabla  1.  (continued) 


OLIVE  GREEN 
SAMPLE 
STANDARD 

THIN  standard 

THIN  YELLOW 
THIN  GREEN 
THIN  RED 
FULL  STANDARD 
FULL  YELLOW 
FULL  GREEN 
full  RED 


107,  NYLON/CTN  <50/50)  POPLINr  5  OZ. 


32.26 

33.59 

33.84 

32.95 
33.52 
31.18 
31.89 
32.04 

31.95 


-3.46 

-3.65 

-3.66 

-3.68 

-3.56 

-3.21 

-3.28 

-3.65 

-3.31 


B  POSITION 

11  37  NOT  ON  BOUNDARY 

11.93  ON  BOUNDARY 

11.95  ON  BOUNDARY 

iJ, -ro  ON  BOUNDARY 

11.74  ON  BOUNDARY 

11. -H  ON  BOUNDARY 

11.73  ON  BOUNDARY 

11.55  ON  BOUNDARY 

11.32  ON  BOUNDARY 


BLUE  15GTR0P.  WL.  10  02. 


SAMPLE 

STANDARD 
THIN  STANDARD 
THIN  BLUE 
THIN  GREEN 
THIN  RED 
FULL  STANDARD 
FULL  BLUE 
FULL  GREEN 
FULL  RED 


15.46 

15.70 

15.61 

15.67 

15.60 

15.70 

15.70 

15.41 

15.58 


B  POSITION 

_4  ON  BOUNDARY 

-4^61  ON  BOUNDARY 

-a  PI  ON  BOUNDARY 

-4 IsO  NOT  ON  BOUNDARY 
-4.60  ON  BOUNDARY 

-4,18  ON  BOUNDARY 

-4.56  ON  BOUNDARY 

-4.24  ON  BOUNDARY 

-4.28  ON  BOUNDARY 


army  green 

SAMPLE 

STANDARD 
THIN  STANDARD 
THIN  BLUE 
THIN  GREEN 
THIN  RED 
FULL  STANDARD 
FULL  BLUE 
FULL  GREEN 
FULL  RED 


3 44,  POLYESTER/UOQL  GABARDINE  9  OZ. 

a  B  POSITION 


27.12 

27.13 
27.30 

27.26 
26.47 

27.27 
26.65 
26.39 
26.44 


-3.52 

-3.39 

-3.42 

-3.60 

-3.47 

-3.13 

-3.50 

-3.68 

-3.19 


ON  BOUNDARY 
NOT  ON  BOUNDARY 
ON  BOUNDARY 
ON  BOUNDARY 
NOT  ON  BOUNDARY 
ON  BOUNDARY 
ON  BOUNDARY 
ON  BOUNDARY 
ON  BOUNDARY 


ARMY  GREEN 

SAMPLE 

S. AN CARD 
THIN  STANDARD 
THIN  BLUE 
THIN  GREEN 
THIN  RED 
FULL  STANDARD 
FULL  BLUE 
FULL  GREEN 
FULL  RED 


344,  POLYESTER/UOOL  TROPICAL,  10  02. 

a  n  POSITION 


27.33 
26.63 
26.37 
27.15 
2/. 10 
25.35 
24.98 
26.07 
25.76 


-3.29 

-3.66 

-3.36 

-3.79 

-3.13 

-3.01 

-3.00 

-3.05 

-2.75 


ON  BOUNDARY 
ON  BOUNDARY 
NOT  ON  BOUNDARY 
ON  BOUNDARY 
ON  BOUNDARY 
ON  BOUNDARY 
ON  BOUNDARY 
ON  BOUNDARY 
ON  BOUNDARY 
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Table  1.  (continued) 


ARMY  GREEN 

44f  UL.  SERGE  15 

OZ. 

SAMPLE 

L 

A 

B 

POSITION 

STANDARD 

26  .  25 

-3.18 

1.11 

ON  BOUNDARY 

THIN  BLUE 

26.10 

-3.53 

1.18 

ON  BOUNDARY 

THIN  GREEN 

26.21 

-3.52 

1.17 

ON  BOUNDARY 

THIN  RED 

26.21 

-3.46 

1.24 

ON  BOUNDARY 

FULL  STANDARD 

26.15 

-3.43 

1.00 

ON  BOUNDARY 

THIN  STANDARD 

26.24 

-3.49 

1  •  22 

ON  BOUNDARY 

FULL  BLUE 

26.47 

-3.46 

.97 

ON  BOUNDARY 

FULL  GREEN 

25.79 

-3.35 

1.09 

ON  BOUNDARY 

FULL  RED 

26.22 

-3.36 

1.13 

NOT  ON  BOUNDARY 

OLIVE  GREEN 

108r  UOOL/NYLONr 

16  OZ. 

SHIRTING 

SAMPLE 

L 

A 

B 

POSITION 

STANDARD 

29.60 

-2.66 

11.77 

ON  BOUNDARY 

THIN  STANDARD 

30.51 

-2.88 

12.35 

ON  BOUNDARY 

THIN  vELLOU 

30.35 

-2.92 

12.58 

ON  BOUNDARY 

THIN  GREEN 

30.52 

-3.41 

12.21 

ON  BOUNDARY 

THIN  RED 

30.19 

-2.74 

12.06 

ON  BOUNDARY 

FULL  STANDARD 

28.52 

-2.75 

11.44 

NOT  ON  BOUNDARY 

FULL  YELLOU 

28.48 

-2.72 

11.64 

ON  BOUNDARY 

FULL  GREEN 

29.28 

-3.16 

11.52 

ON  BOUNDARY 

FULL  RED 

27.81 

-2.41 

10.70 

ON  BOUNDARY 

OLIVE  GREEN 

107.  CTN 

.  SATEEN 

!»  8.2  OZ 

• 

SAMPLE 

L 

A 

B 

POSITION 

STANDARD 

31.30 

-3.46 

10.31 

DN  BOUNDARY 

THIN  STANDARD 

31.78 

••4.11 

11.10 

ON  BOUNDARY 

THIN  YELLOU 

31.78 

-3.86 

11.58 

ON  BOUNDARY 

THIN  RED 

31.60 

-3.58 

11.58 

ON  BOUNDARY 

THIN  GREEN 

30.64 

-4.13 

10.84 

ON  BOUNDARY 

FULL  STANDARD 

30.84 

-3.56 

10.74 

NOT  ON  BOUNDARY 

FULL  YELLOU 

30.62 

-3.14 

10.84 

ON  BOUNDARY 

FULL  GREEN 

•  30.35 

->■73 

10.70 

ON  BOUNDARY 

FULL  RED 

31.53 

-3.04 

10.88 

ON  BOUNDARY 

BLUE  151 »  TROP.  WL.r 

1  :•  t  7 

SAMPLE 

L 

A 

B 

POSITION 

STANDARD 

23.48 

4.48 

-17.33 

NOT  ON  BOUNDARY 

THIN  STANDARD 

23.85 

4.68 

-17.47 

ON  BOUNDARY 

THIN  BLUE 

23.54 

4.39 

-17.10 

ON  BOUNDARY 

THIN  GREEN 

24.44 

4.65 

-18.13 

ON  BOUNDARY 

THIN  RED 

24.18 

4.61 

-17.72 

ON  BOUNDARY 

FULL  STANDARD 

23. C7 

4.40 

-17.07 

NOT  ON  BOUNDARY 

FULL  BLUE 

22.84 

4.34 

-16.72 

FULL  GREEN 

23 . 26 

4.43 

-17.52 

ON  .  ..  RY 

FULL  RED 

23.12 

4.62 

-17.19 

ON  BOUNDARY 

11 
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Figure  2 
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L*  vs.  a*  plot  of  Standard  and  limit  samples  for  Olive  Drab  7 
untreated  cotton  duck,  8.25  oz. 
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color  space,  we  wrote  a  computer  program  that  determined  the  minimum  nor.- 
tilting  ellipsoid  that  could  be  constructed  around  these  samples.  The 
non-tilting  feature  meant  that  the  ellipsoid  had  a  horizontal  major  axis, 
a  horizontal  minor  axis,  both  in  the  chromaticity  plane,  and  a  vertical 
axis;  the  two  horizontal  axes  lay  entirely  in  the  a*-  b*  plane  and  the 
vertical  axis  was  perpendicular  to  the  other  two.  We  propose  this 
non-tilting  restriction  in  the  spirit  of  the  observations  of  Schultze 
and  Gall  (l)  and  of  MacAdam  (2),  as  applied  to  perceptibility  ellipsoids. 

Each  ellipsoid  was  therefore  defined  by  seven  parameters;  the  three 
position  coordinates  of  the  center,  the  length  of  the  three  semiaxes,  and 
the  angle  tnat  the  major  axis  in  the  chromaticity  plane  makes  with  the  a* 
axis.  The  computer  program  varied  these  seven  parameters  one  at  a  time 
in  a  systematic  fashion,  each  time  determining  whether  a  smaller  ellipsoid 
could  be  constructed  that  contained  all  the  points  for  the  standard  and 
limit  samples.  The  program  stopped  when  no  parameter  ”ariation  was  possible 
that  would  shrink  the  ellipsoid  without  excluding  one  of  the  points. 

Appendix  B  describes  the  mathematics  used  and  presents  the  computer  program, 
called  ELPSOID.  The  results  from  this  program  are  given  in  Table  2.  In 
this  table  the  12  colors  are  presented  in  the  same  order  as  they  are  in 
Table  1. 

2.  Non-tilting  ellipsoid,  axis  angle  equal  to  hue  angle.  We  wrote 
another  minimum  ellipsoid  program,  similar  to  the  first  except  that  this 
time  we  had  only  six  adjustable  parameters  instead  of  seven.  The  axis 
angle  was  set  equal  to  the  hue  angle;  the  latter  was  the  angle  between 
the  line  connecting  the  center  of  the  ellipsoid  with  the  L*  axis  and  the 
a*  axis.  Thus,  one  of  the  horizontal  ax«s  of  the  ellipsoid  always  coincides 
with  the  isohue  direction,  the  other  horizontal  axis  with  the  isochroma 
direction.  Figure  3  serves  as  an  illustration,  in  two-dimensional  space, 
of  the  difference  between  the  two  programs.  The  color  points  shown  were 
plotted  arbitrarily.  The  two  ellipses  drawn  on  the  figure,  however,  were 
calculated  by  a  variation  of  the  two  programs  that  constructs  ellipses  instead 
of  ellipsoids.  The  dashed  ellipse  is  the  smallest  ellipse  of  any  orientation 
that  can  be  drawn  surrounding  all  the  points  whose  major  axis  coincides  with 
an  isohue  line. 

The  reason  for  this  second  minimum  ellipsoid  program  was  that  we  wished 
to  take  advantage  of  the  commercial  experience  which  nas  been  built  up  in  the 


1.  W.  Schultze  and  L.  Gall,  ’’Experimentelle  Uberprufung  Mehrerer 
Farbabstandsformelu  bez'^glich  du  Helligkeits  -  und  Sattigungsdifferenzen 
bei  gesattigten  Farben,"  Farbe  18,  131-148  (1969). 

2.  i).  L.  MacAdam,  in  Color  Metrics,  edited  by  ,T.  J.  Vos,  L.F.C,  Friele  and 
P.L.  Walraven  (AIC/Holland,  Soesterburg,  1972)  pp  160-170. 


Table  2.  Minimum  Ellipsoid  Characteristics 


COLOR  NO. 


10.84  12.08  15.73 


1  2  3 

Unconstrained  Azimuthal  Angle 
Ellipsoid  , 

Center:  L  34.44  G8.91  54.89 
a*  -2.42  .77  2.8 

b*  10.84  12.98  15.7 

Semiaxis 

Lengths: 

Maj.Horiz.  1.58  .99  .4 

Min.Horiz.  .26  .65  .4 

Vertical  2.56  2.62  1.2 

<  botween 
Horiz.axis  k 

a  axis, deg.  -48.2  -62.0  73.6 
Vol.  4.40  7.07  1.0 

Constrained  Azimuthal  Angle 


Ellipsoid  , 
Center:  L 

* 

a 

* 

b 

Semiaxis 

Lengths: 

Maj.Horiz. 

Min.Horiz. 

Vertical 

<  between 
Horiz.axlB  k 
a*  axis. deg. 

Vol. 


34.48  68.79  54.81 


10.86  13.18  15.77 


-77.4 

19.04 


4 

5 

6 

7 

8 

9 

10 

U 

12 

26.25 

32.51 

15.60 

26.84 

26.27 

26.24 

29.11 

31.23 

23.64 

-3.33 

-3.43 

1.13 

-3.41 

-3.27 

-3.36 

-2.86 

-3.66 

4.49 

1.18 

11.58 

-  .44 

1.31 

1.01 

1.10 

11.74 

1U.97 

-17.42 

.53 

.55 

.48 

.39 

.61 

.19 

1.35 

.78 

1.17 

.27 

.36 

.17 

.19 

.50 

.19 

.47 

.60 

.19 

.55 

2.15 

.22 

.63 

1.78 

.45 

2.40 

.99 

1.05 

68.8 

-58.7 

-76.4 

1.8 

-23.3 

a 

-65.8 

-55.7 

-80.4 

.33 

1.78 

.08 

.19 

2.28 

.07 

6.36 

1.95 

.98 

26.22 

32.44 

15.60 

26.89 

26.40 

26.21 

29.47 

31.16 

23.65 

-3.35 

-3.50 

1.15 

-3.43 

-3.23 

-3.37 

-2.85 

-3.63 

4.51 

1.26 

.1.47 

-4.49 

1.33 

1.03 

1.07 

11.41 

10.96 

-17.36 

.2lt 

.66 

.44 

.45 

.63 

.22 

1.54 

.91 

1.06 

•  60b 

.31 

.15 

.22 

.43 

.17 

.52 

.62 

.16 

.63 

2.21 

.24 

.69 

2.04 

.43 

2.37 

.90 

1.37 

•20.6 

-73.0 

-75.6 

-21.1 

-17.7 

-17.7 

-76.0 

-71.7 

-75.4 

.34 

1.86 

.07 

.28 

2.32 

.07 

7.94 

2.14 

.99 

'  Two  horizontal  axes  equal;  therefore  angle  cannot  be  specified. 

b  Iso(hue-lightne8s)  axis  (considered  by  computer  program  to  be  major  horizontal  axis)  is  shorter  than 
iso(chroma-llghtness)  axis. 
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HUE  ANGLE 

_ I _ i 

-1.5  -1.0 


the  two  minimmi  ellipsoid. 


textile  industry  relating  to  acceptance  of  production  lots.  This  experience 
has  indicated  that  acceptability  tolerances  in  a  given  uniform  color  space 
should  be  so  arranged  as  to  allow  the  greatest  leeway  in  the  lightness 
direction,  intermediate  latitude  in  the  chroma  direction,  and  the  most 
stringent  tolerances  in  the  hue  direction  (3).  We  can  thus  picture  an 
acceptability  ellipsoid  with  the  largest  axis  pointing  in  the  L*  direction, 
the  intermediate  axis  in  the  chromaticity  plane  along  an  isohue  line,  and 
the  shortest  axis  in  the  chromaticity  plane  perpendicular  to  the  isohue 
line.  The  ratios  between  these  three  tolerances  have  typically  been  found 
to  range  between  3:2:1  and  4:2:1,  in  CIELAB  space. 

This  computer  program  was  called  ELPS0D2,  and  is  given  in  Appendix  B. 

Results  of  this  program  are  presented  in  Table  2. 

C.  Results  of  ellipsoid  tests;  hue- chroma- value  axis  ratios. 

Comparison  of  the  axis  angles  given  by  the  two  minimum  ellipsoid 
programs  is  interesting  and  is  shown  in  Table  3,  which  was  abstracted  from 
the  results  in  Table  2.  We  see  a  remarkable  agreement  between  the  two 
programs  with  respect  to  the  axis  angle  of  the  ellipsoid.  In  six  out  of 
eleven  cases  the  angles  agreed  to  10  or  less;  in  eight  out  of  eleven 
cases  the  agreement  is  to  20  or  less.  Since  in  one  of  the  programs, 
one  horizontal  axis  was  constrained  to  the  isohue  direction,  whereas  in 
the  other  program  the  axes  were  unconstrained,  we  conclude  that  the 
acceptability  ellipsoids,  at  least  insofar  as  they  are  defined  by  the 
standard  and  limit  samples,  naturally  tend  to  orient  themselves  in  the 
isohue  and  isochroraa  directions. 

Table  4,  calculated  from  the  figures  in  Table  2,  shows  the  ratios  of 
the  three  principal  axes  for  the  constrained  ellipsoids,  as  well  as  the 
ellipsoid  volumes.  We  see  that  in  nine  out  of  the  twelve  cases  we  have  the 
expected  size  relationship  of  the  largest  axis  being  in  the  lightness 
direction,  the  intermediate  axis  in  the  chroma  direction,  and  the  smallest 
axis  in  the  hue  direction.  In  one  of  the  other  three  cases  the  lightness 
and  chroma  axes  were  equal;  in  the  second  of  the  three  exceptions  the 
chroma  and  hue  axes  were  reversed  from  the  expected  order,  and  in  the 
third  case  the  lightness  and  chroma  axes  were  reversed.  These  results 
seem  to  bear  out  the  commercial  experience  for  textile  industry  color 
acceptances  mentioned  above. 

III.  TESTS  ON  MULTIPLE  SAMPLE  COLLECTIONS 
A.  Selection  of  sample  pairs 

As  a  gtneral  policy,  we  thought  that  acceptability  standards  should  be 
based  on  the  visual  judgment  of  trained  and  experienced  observers.  In  view  of 
the  ellipsoid  orientation  indicated  by  the  work  just  described,  we  thought  that 
it  would  be  appropriate  to  establish  six  tolerances:  two  in  the  plus  ana  minus 

3.  "The  Hue  Angle  Sort  Program,"  pamphlet  issued  by  Hunter  Associates 
Laboratory,  Inc. 
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TABLE  3.  ANGLE  OF  PRINCIPAL  HORIZONTAL  AXIS  OF 
ELLIPSOID  WITH  a*  COORDINATE  AXIS,  IN 
DEGREES. 


Olive  Drab  7,  Cotton  Duck 
Olive  Green  107,  Nylon/cotton 
Olive  Green  108,  Wool/nylon 
Olive  Green  107,  Cotton  sateen 
Army  Green  44,  Wool  Gabardine 
Army  Green  344,  Polyester/Wool 
Gabardine 

Army  Green  344,  Tropical  Wool 

Army  Green  44,  Wool  Serge 

Tan  46,  Cotton  Poplin 

Tan  M-l,  Polyester/Wool  Tropical 

Blue  150,  Tropical 

Blue  151,  Tropical 


Angle 

unconstrained 

Angle 

constrained 
along  isohue 

1  ine 

Corrected 

Difference 

-48 

-77 

29 

-59 

-73 

14 

-$6 

-76 

10 

-56 

-72 

16 

69 

-21a 

0 

2 

-21 

23 

-23 

-18 

5 

b 

-18 

-62 

86 

32 

74 

80 

6 

-76 

-76 

0 

-80 

-75 

5 

a  Axis  in  isohue  direction  smaller  than  axis  in  isochroma  direction;  angle 
must  therefore  be  increased  by  90°  before  comparison  with  value  in  left 
column. 

b.  Two  horizontal  axes  equal;  therefore  angle  cannot  be  specified. 

TABLE  AXIS  RATIOS  AND  ELLIPSOID  VOLUMES  FOR 

ELLIPSOIDS  POINTING  IN  ISOHUE  DIRECTION 


Lightness. Chroma s  Hue  Ellipsoid 

Volume 


Olive  Drab  7,  Cotton  Duck 
Olive  Green  107,  Nylon/ Cotton 
Olive  Green  108,  Wool/Nylon 
Olive  Green  107,  Cotton  Sateen 
Army  Green  44,  Wool  Gabardine 
Army  Green  344,  Polyester/Wool 
Army  Green  344,  Tropical 
Army  Green  44,  Wool  Serge 
Tan  46,  Cotton  Poplin 
Tan  M-l,  Polyester/Wool  Tropical 
Blue  150,  Tropical 
Blue  151,  Tropical 


3:37  s  1.66  :  1.00  19.04 

7.23  s  2.16  S  I. 00  1.86 

4.59  J  2.98  s  1.00  7.94 

1.45  :  1.45  s  1.00  2.14 

1.05  :  0.36  s  1.00  0.34 

3.23  :  2.10  :  1.00  0.28 

4.72  :  1.45  :  1.00  2.32 

2.56  s  1.29  :  1.00  0.07 

3.44  :  1.72  :  1.00  9.96 

3.00  :  1.10  :  1.00  1.12 

1.62  I  3.01  s  1.00  0.07 

8.41  3  6.53  8  1.00  0.99 
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chroma  directions,  two  in  the  plus  and  minus  lightness  directions,  and  two  in 
the  plus  and  minus  hue  directions.  The  tolerances  would  be  set  up  on  the  basis 
of  the  responses  of  the  inspectors  to  samples  shown  to  them.  We  would  then 
establish  an  acceptability  ellipsoid  that  incorporates  these  tolerances  and 
that  has  its  three  axes  oriented  in  the  hue,  chroma,  and  lightness  directions. 

The  problem  then  arises  as  to  how  to  obtain  samples  that  differ  from  the 
standard  only  in  hue,  chroma,  or  lightness  and  in  varying  amounts  for  each. 

We  rejected  a  preliminary  plan  involving  the  dyeing  of  samples  of  various 
degrees  of  color  difference  in  each  of  the  six  directions  when  we  realized 
how  difficult  it  would  be  to  prepare  such  dyeings.  In  order  to  obtain 
samples  having  the  precise  location  in  color  space  that  we  want,  wt  would 
have  to  have  at  least  two  and  probably  more  trial  dyeings  for  each  sample. 

Such  a  procedure  would  take  an  inordinate  amount  of  time. 

A  much  better  plan  seemed  to  be  to  try  to  find  the  desired  samples  among 
previous  submissions  of  samples  from  textile  mills.  For  each  shade  that  is 
accively  used  by  the  Army,  textile  mills  are  continually  submitting  samples 
that  are  either  accepted  or  rejected  by  the  inspectors.  Fortunately,  a  set  of 
some  200  to  300  samples  was  available  for  each  of  three  shades  of  interest  to 
the  Arrays  olive  green,  tan,  and  dark  blue.  These  samples  are  located  in  a 
cloud  of  points  close  to  the  standard,  since  they  are  intended  to  match  the 
standard.  It  should  be  possible  to  find,  among  all  these  samples,  test  cases 
that  could  be  shown  to  the  observers  and  that  would  serve  to  define  the  re¬ 
quired  tolerances. 

We  wished  to  find  samples  representing  four  color  differences  of  increasing 
amount  in  each  of  the  six  directions  mentioned  above.  These  color  differences 
would  ideally  extend  from  completely  acceptable  (for  the  smallest  color  difference) 
to  completely  unacceptable  (for  the  largest  color  difference).  The  two  intermediate 
color  differences  would  presumably  be  close  to  the  tolerance  limits.  The  procedure 
would  be  to  show  each  of  the  color  differences  to  an  observer  ten  times  (randomly 
mixed  with  other  samples  so  that  the  inspector  does  not  realize  that  he  is  looking 
at  the  same  sample  repeatedly).  Each  time  that  the  observer  is  shown  the  color 
difference,  he  is  asked  whether  he  would  accept  or  reject  such  a  difference  if  it 
were  an  incoming  shipment  against  standard.  The  experimental  regime  is  known 
as  the  method  of  constant  stimuli. 

To  explain  the  procedure,  let  us  assume  that  the  four  color  differences 
chosen  are  0.5,  1.0,  1.5,  and  2.0  CIELAB  units.  Let  us  observe  that  the  ob¬ 
server  accepts  the  0.5  color  difference  ten  times  out  of  ten;  the  1.0  color 
difference  seven  times  out  of  ten;  the  1.5  color  difference  three  times  out 
of  ten.  It  can  then  easily  be  seen  that  the  507.  acceptance  limit  should  be 
close  to  1.25  color  difference  units.  In  fact,  in  the  more  general  case, 
the  method  of  logistic  functions  (see  below)  is  used  to  determine  the  507. 
acceptability  limit. 


The  question  then  arises  as  to  how  to  select  these  color  differences  out  of 
the  cloud  of  200  to  300  samples.  We  realize  that  it  is  not  necessary  to  have  only 
a  single  standard  against  which  all  the  judgments  would  be  nade.  All  that  is 
necessary  is  a  selection  of  sample  pairs,  of  which  one  ould  play  the  role  of 
standard  and  the  other  would  play  the  role  of  sample  submitted  for  acceptance. 
Since  we  have  six  directions  in  color  space  with  four  degrees  of  color  differ¬ 
ence  for  each  direction,  we  need  24  sample  pairs.  These  pairs  were  chosen  in 
the  following  ways 

We  first  measured  all  the  samples  for  each  of  the  three  colors  on  the 
Diano-Hardy  spectrophotometer,  and  obtained  the  CIELAB  coordinates  (D75,  CIE 
1931).  We  then  created  a  computer  file  with  an  entry  for  each  of  the  200  to  300 
samples  for  each  color;  the  file  contained  the  numerical  sample  designations 
and  the  three  C1ELAB  coordinates  L*a*b*.  We  rearranged  this  file  so  that  it 
was  in  increasing  numerical  order  for  the  sample  designations.  This  rear¬ 
rangement  was  carried  out  by  Program  REARRAN,  given  in  Appendix  C. 

We  then  wrote  a  computer  program  which  we  called  Program  SIXWAYS.  For 
each  of  the  three  collections  of  samples,  the  program  took  each  individual 
sample,  one  at  a  time,  and  constructed  three  lines  in  CIELAB  color  space 
passing  through  the  sample.  The  first  line  was  an  iso  (hue-lightness)  line, 
which  means  that  it  was  a  line  in  the  a*  -  b-'  plane  that  joined  the  origin 
with  the  point  in  question.  Since  the  line  passed  through  the  point  for  the 
sample,  it  extended  from  the  sample  in  two  directions:  the  direction  of  in¬ 
creasing  chroma  and  the  direction  of  decreasing  chroma.  The  second  line  was 
an  iso  (chroma- lightness)  line,  or  a  line  in  the  a*  -  b*  plane  that  was  per¬ 
pendicular  to  the  first.  This  line  also  passed  through  the  sample  and  ex¬ 
tended  from  the  sample  in  two  opposite  hue  directions.  The  third  line  was 
an  iso  (hue-chroma)  line,  or  a  vertical  line  through  the  sample  perpendicular 
to  the  a*  -  b*  plane  and  extending  from  the  sample  in  a  plus  lightness  and  a 
minus  lightness  direction. 

The  computer  program  thus  contained  an  outer  loop  that  was  traversed  n 
times,  where  n  is  the  number  of  samples  in  the  collection.  For  each  of  these 
samples,  the  program  constructed  the  six  lines  in  CIELAB  color  space  just  de¬ 
scribed.  The  program  then  went  through  an  inner  loop,  with  n  -  1  passages,  during 
which  it  again  considered  every  sample  in  the  collection  and  determined  whether  it 
lay  on  one  of  the  six  lines  that  was  constructed  (the  program  actually  determined 
whether  the  sample  deviated  from  these  lines  by  no  more  than  0.1  CIELAB  unit). 

If  the  sample  did  so  lie-  the  program  calculated  the  color  difference  between 
the  two  samples.  We  thus  had,  fer  each  of  the  samples  in  the  collection  con¬ 
sidered  temporarily  as  a  standard,  a  subset  of  samples  extending  from  the 
standard  in  each  of  the  six  directions,  together  with  the  associated  color 
differences  between  each  sample  in  the  subset  and  the  temporary  standard. 

It  was  therefore  relatively  easy  to  select  pairs  of  samples  for  which  the 
"sample"  deviated  from  the  "standard"  in  each  of  six  directions  an<*  with  four 
different  color  differences  in  each  of  the  directions.  Notice  that  there  was  no 
single  fixed  standard  chosen;  the  standard  varies  with  each  pair.  Notice  also 
that  the  sizes  of  the  color  differences  available  were  determined  and  limited 
by  the  particular  samples  in  the  collection. 
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Appendix  D  shows  the  mathematics  that  was  used  to  determine  how  far  a 
given  sample  deviated  from  the  six  iso  lines  radiating  from  another  sample. 

The  same  appendix  presents  the  listing  of  programs  SIXWAYS.  Also  included 
in  the  appendix  is  a  typical  page  from  the  output  of  program  SIXWAYS. 

Figure  4  shows  an  a*  -  b*  plot  of  the  samples  making  up  the  four  color 
differences  of  increasing  chroma  and  the  four  color  differences  of  increasing 
hue  for  the  olive  green  collection.  The  points  represent  actual  samples  in 
the  collection;  the  solid  lines  join  pairs  of  points  for  the  various  color 
differences.  The  dashed  lines  are  the  true  constant  chroma  lines  and 
constant  hue  lines,  and  it  can  be  seen  that  in  no  case  does  a  sample 
deviate  from  these  lines  by  more  than  0.1  CIELAB  unit. 

We  were  able  to  physically  reduce  the  number  of  p4lrs  of  samples 
chosen  by  almost  half,  because  in  all  but  a  few  cases,  the  same  pair  of 
samples  was  used  for  a  minus  difference  as  for  a  plus  difference.  For 
example,  let  us  suppose  that  Samples  A  and  B  have  the  same  hue  and  light¬ 
ness  but  differ  by  1.0  CIELAB  unit  in  chroma.  In  order  to  determine 
whether  there  is  a  difference  in  acceptability  in  the  plus  chroma  or  the 
minus  chroma  direction,  the  pair  will  be  shown  to  the  inspector  in  two 
different  ways;  in  the  first  way.  Sample  A  is  the  standard  and  Sample  B 
is  the  sample  to  be  judged  against  the  standard;  in  the  second  way.  Sample 
B  is  the  standard  and  Sample  A  is  the  sample.  In  fact,  if  Sample  B  would 
appear  under  the  Sample  A  heading  in  the  printout  from  progrf.u  SIXWAYS  as 
being  more  chromatic  than  Sample  A,  then  Sample  A  would  appear  under  the 
Sample  B  heading  as  being  less  chromatic. 

B.  Presentation  of  Sample  Pairs 

As  a  result  of  work  described  in  the  preceding  section,  we  had  72 
sample  pairs  to  show  to  the  observers  (three  shades,  six  directions  in  color 
space  for  each  shade,  that  is  4-  chroma,  +  lightness,  and  four  degrees  of 
color  difference  for  each  direction).  We  wished  to  show  each  sample  pair 
to  each  observer  ten  times,  making  a  total  of  720  Individual  observations 
by  each  inspector.  In  order  not  to  confuse  the  observers  more  than  was 
absolutely  necessary,  we  did  one  shade  at  a  time,  involving  240  judgments 
per  shade.  We  wrote  a  computer  program,  entitled  Program  FUNFER,  that  used  a 
random  number  generator  to  randomize  the  order  of  presentation  of  these  240 
trials  (see  Appendix  E).  The  observers  were  not  told  that  there  were  replicate 
sample  presentations,  although  they  probably  guessed  that  from  the  fact  that 
there  were  many  fewer  sample  pairs  chan  there  were  presentation.  In  order  to 
prevent  the  inspector  being  influenced  by  the  location  of  the  sample  pair  as 
he  saw  it  being  removed,  we  had  a  "lazy  susan"  sample  holder  constructed,  on 
which  the  samples  were  arranged  circularly  and  accessed  by  rotating  a  tray 
In  this  way,  distinctive  sample  positions  were  abolished. 

We  administered  the  test  to  six  observers,  three  from  the  government 
service  and  three  from  industry.  In  addition,  a  seventh  observer,  previously 
from  industry  but  now  in  the  government  service,  looked  at  only  one  of  the 
three  shades— the  dark  blue.  The  observer  was  presented  with  a  sample  pair 
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Figure  4.  Sample  pairs  for  Olive  Green  selected  for  observers’  evaluation 
plotted  on  a*  -  b*  diagram. 
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and  was  told  which  was  the  standard  and  which  was  the  sample  to  be  compared  with 
standard.  He  or  she  was  told  to  either  accept  or  reject  the  sample  against 
standard;  no  intermediate  judgments  were  allowed.  The  observers*  judgments 
were  tabulated,  and  the  percentage  of  passes  for  each  pair  v;as  noted. 

Tables  5,  6  and  7  show  the  raw  results  of  the  observers  for  all  the 
samples.  In  the  leftmost  column,  the  CIELAB  total  color  difference  for  eaoh 
pair  is  given.  In  the  columns  to  the  right  the  percent  pass  judgments  for 
each  pair  is  shown.  A,  B  and  C  were  the  three  observers  from  the  government 
service;  D,  E,  and  F  were  the  three  observers  from  industry.  G  was  the  ob¬ 
server  who  was  previously  in  in  industry  and  is  now  in  the  government  service, 
who  looked  at  only  the  dark  blue  series  of  samples. 


We  see  that  in  most  cases  the  percentage  of  samples  passed  decreases 
with  increasing  color  difference.  Occasionally  there  are  reversals,  as  in¬ 
dicated  by  a  larger  color  difference  giving  a  greater  percentage  of  acceptances. 
In  some  cases  these  reversals  are  severe  enough  to  invalidate  that  particular 
set  of  judgments,  as  will  be  seen  later.  For  the  most  part,  however,  the  re¬ 
sults  of  these  tests  are  surprisingly  good. 

C.  Scoring  by  Logistic  Function 

The  raw  data  in  Tables  5,  6  and  7  take  the  form  of  percent 
acceptance  as  a  function  of  color  difference.  We  now  have  to  calculate  from 
these  data,  on  a  valid  statistical  basis,  the  color  difference  that  would 
correspond  to  a  50%  acceptance  rate.  We  used  the  method  of  logistic  functions, 
explained  in  Appendix  F.  This  method  not  only  gives  the  50%  point  but  also 
gives  some  idea  of  the  precision  of  this  determination  in  the  form  of  a 
standard  deviation. 

In  Tables  5f  6  and  7  '  we  see  next  to  each  group  of  four  acceptance 

percentages  the  calculated  color  difference  corresponding  to  a  50%  acceptance 
rate,  together  with  the  standard  deviation  of  this  color  difference.  This 
standard  deviation  is  really  a  measure  of  the  self-consistency  of  the  acceptance 
judgments.  If  there  are  serious  reversals  in  the  judgments,  the  standard  de¬ 
viation  may  be  quite  large  in  comparison  with  the  value  itself  (in  some  cases 
larger). 

We  decided,  quite  pragmatically,  to  accept  only  those  values  that  were  at 
least  about  four  times  as  large  as  their  standard  deviations.  In  addition,  in 
view  of  the  fact  that  observer  F  accepted  18  out  of  the  24  tan  pairs  and  20 
out  of  the  24  dark  blue  pairs  on  a  100%  basis,  we  considered  him  unfamiliar 
with  colors  and  used  his  results  only  for  the  olive  green  pairs. 

D.  Results  of  Tests 

Table  8  summarizes  the  tolerance  results  that  were  accepted.  The 
standard  deviations  shown  refer  to  the  values  given  in  this  table  and  do  not 
involve  the  standard  deviations  shown  in  Tables  5,  6  and  7.  The  following 
observations  can  be  made* 
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Table  5.  Raw  Results  for  Observers  A,  B,  and  C,  Together  with  50% 
Acceptance  Limits 

color:  olive  green 


DELTA 

c 

A 

B 

C 

n 

PCT . 

BOUNDARY 

PCT .  BOUNDARY 

PCT.  BOUNDARY 

PASS 

t  STD.  DEV. 

PASS  1  STD.  DEV. 

PASS  1  STD.  DEV. 

PLUS  CHROMA 

.40 

100 

100 

100 

.80 

40 

50 

70 

1.13 

20  .78  ± 

.11  80 

1.31  ±  .46  90 

4.02  ±  7.01 

1.52 

20 

40 

70 

MINUS  ChROMA 


.40 

100 

100 

100 

.80 

100 

50 

80 

i.13 

100 

0.00  i  0.00  70 

1.05  ±  .19  80 

1.47  ±  .38 

1 .52 

100 

20 

40 

PLUS  HUE 

.40  100  70  100 

.77  70  0  40 

1.17  10  .83  ±  .08  0  .46  ±  .07  0  .71  ±  .08 

1.32  0  0  0 


MINUS  HUE 

.40  40  90  90 

.77  20  0  0 

1.17  0  .34  ±  .16  0  .53  ±  .07  0  .53  ±  .07 

1.45  0  0  0 


PLUS  LIGHT 


.60 

too 

90 

100 

1.17 

70 

100 

100 

1.82 

30 

1.76  ±  .38  10 

1.64  ±  .3,4  70 

6.70  ±13.73 

2.29 

50 

40 

90 

MINUS  LIGHT 


.60 

100 

100 

100 

1.17 

100 

ICO 

100 

1.82 

100 

4.46  ±  3.11  10 

1.78  ±  .29  100 

2.29  ±  .15 

2.29 

90 

40 

50 
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Table  5.  (continued) 


COLORS 


DELTA  A 

E 

PCT.  BOUNDARY  PCT. 
PASS  l  STD.  DEV.  PASS 


PLUS  CHROMA 


.50 

80 

100 

.92 

80 

40 

1.33 

10 

.95  4  .15  0 

1.92 

10 

0 

HINUS  CHROMA 

.50 

90 

100 

.92 

20 

20 

1.33 

0  .73  4 

.10  10 

1.92 

0 

0 

PLUS  HUE 

.25  40  70 

.47  40  60 

.69  0  .21  4  .17  0 

1.01  0  0 


MINUS  HUE 


.25  50  70 

.48  10  0 

.69  0  .25  4  .07  0 

1.01  0  0 


PLUS  LIGHT 


.66 

60 

100 

1.32 

30 

50 

1.95 

20 

.84  4  .25  0 

2.65 

0 

0 

MINUS 

LIGHT 

.66 

100 

100 

1.32 

50 

20 

1.95 

30 

1.41  4 

.17  70 

2.65 

10 

0 

TAN 


B  C 

BOUNDARY  PCT .  BOUNDARY 
t  STD.  DEV.  PASS  l  8TD.  DEV. 


80 

80 

.85  6  .09  20  .98  4  .15 

10 


100 

30 

.?7  4  .13  0  .80  4  .10 

0 


60 

30 

.43  4  .10  0  .30  4  .08 

0 


60 

7C 

.29  4  .04  10  .42  4  .15 

0 


100 

20 

1.27  4  .14  10  1.09  *  .18 

0 


100 

30 

1.50  4  .38  10  1.15  .  .16 

0 
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rable  5.  (continued) 


COLORS  DARK  BLUE 


DELTA 

E 

PCT  • 
PASS 

A 

BOUNDARY 
t  STD.  DEV. 

PCT. 

PASS 

B 

BOUNDARY 

4  STD.  DEV. 

PCT. 

PASS 

G 

BOUNDARY 
l  STD.  DEV. 

PLUS  CHROMA 

.36 

.69 

1.07 

1.42 

70 

20 

0 

20 

.47  ±  .14 

90 

80 

0 

70 

1.77  ±  1.17 

80 

80 

10 

40 

.93  ±  .23 

MINUS  CHROMA 

.36 

.69 

1.07 

1.42 

80 

30 

10 

30 

.59  ±  .15 

100 

90 

100 

80 

6.64  ±14.74 

40 

80 

10 

10 

.42  ±  .22 

PLUS  :*ue 

.18 

.38 

.57 

.75 

70 

80 

60 

50 

1.12  ±  1.29 

100 

100 

100 

90 

1.67  ±  1.48 

100 

90 

80 

0 

.69  ±  .20 

MINUS 

HUE 

.18 

.38 

.57 

.75 

70 

50 

100 

60 

S3. 62  ±*1.01 

100 

90 

100 

10 

.61  ±  .14 

80 

90 

100 

0 

.90  ±  .68 

PLUS  LIGHT 

.50 

1.00 

1.50 

2.01 

50 

30 

0 

0 

.53  ±  .20 

100 

100 

10 

0 

1.29  ±  .12 

60 

10 

0 

0 

.56  ±  .13 

MINUS 

light 

.50 

1.00 

1.50 

2.01 

90 

70 

10 

0 

1.04  ±  .13 

100 

80 

:  20 
0 

1.17  ±  .13 

70 

10 

0 

0 

.62  ±  .12 
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Table 


6,  Raw  Results  for  Observers  D,  E,  end  F,  Together  with  SOT 
Acceptance  Limits 


color:  olive  green 


DELTA 

E 

PCT. 

PASS 


PLUS  CHROMA 


.40 

.80 

1.13 

1 .52 

100 

90 

SO  1.38  k 

50 

.32 

100 

90 

40 

30 

1.13  k 

100 
100 
.17  100 

30 

1.42  k  .10 

MINUS  CHROMA 

.40 

.80 

1.13 

1.52 

100 

40 

100  .92  k 

20 

.18 

100 

60 

20 

10 

.85  k 

100 
100 
.09  70 

10 

1.21  k  .11 

PLUS  HUE 

.40 

.77 

1.17 

1.32 

80 

0 

0  .49  k 

0 

.07 

100 

0 

0 

0 

.55  k 

100 

70 

.08  10 

0 

.83  k  .08 

MINUS  HUE 

.40 

.77 

1.17 

1.45 

100 

0 

0  .55  k 

0 

.08 

100 

0 

0 

0 

.55  k 

100 

10 

.08  0 

0 

.60  k  .10 

PLUS  LIGHT 

•  60 
1.17 
1.82 
2.29 

100 

100 

0  1.46  k 

0 

.12 

100 

90 

30 

0 

1.50  k 

100 
100 
.18  40 

90 

3.15  k  3.52 

MINUS  LIGHT 

.60 

1.17 

1.82 

2.29 

10  0 

100 

0  1.58  ± 

10 

.17 

100 

90 

20 

0 

1.42  k 

100 

100 

.16  60 
90 

4.97  k  8.13 

BOUNDARY 
i  aTD.  DEV. 


PCT. 

PASS 


boundary 

*  am.  DFU . 


PCT. 

PASS 


BOUNDARY 
t  STD.  DEV. 
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Table  6.  (continued) 


color:  tan 


DELTA  u 

E 

PCT.  BOUNDARY  PCT. 
PASS  *  STD.  DEV.  PASS 


E  F 

BOUNDARY  PCT.  BOUNDARY 
S  STD.  DEV.  PASS  t  STD.  DEV. 


PLUS  CHROHA 

.50  100 

.92  0 

1.33  100  1.04  ± 

1.92  0 


100 

90 

26  100*667.20  ±*3.07 

100 


100 

100 

80*378.07  ±*3.62 
100 


MINUS 

CHROMA 

.30 

100 

100 

.92 

0 

100 

1.33 

70 

1.22  ± 

.28  JOO 

1.92 

0 

SO 

100 

100 

1.92  ±  .20  100  5.60  ±  7.94 

90 


PLUS  HUE 


.23 

.47 

.69 

1.01 

100 

90 

0 

0 

.54  ± 

100 

90 

.08  20 

0 

.57  ± 

100 

100 

.07  100*076.48 

100 

±*8.40 

MINUS 

HUE 

.25 

.48 

.69 

1.01 

100 

IOC 

10 

0 

.59  ± 

100 

90 

.04  10 

0 

.55  ± 

10 

100 

,07  20  .65 

100 

±  .18 

PLUS  LIGHT 

.66  100 

1.32  20 

1.95  0  1.07  ± 

2.65  0 


100 

100 

16  100  2.55  ±  .20 

40 


100 

100 

100  2.47  ±  .17 

30 


MINUS  LIGHT 

.66  100 

1.32  0 

1.95  0  .93  ± 

2.65  0 


100 

90 

14  100  2.64  ±  .71 

40 


100 

100 

100  6.47  ±  6.89 

90 
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Table  6.  (continued) 


color:  dark  blue 


DELTA 

D 

E 

E 

PCT .  BOUNDARY 

PASS  t  STD.  DEV. 

PCT. 

PASS 

BOUNDARY 
t  STD.  DEV. 

F 

PCT.  BOUNDARY 
PASS  t  STD.  DEV . 


PLUS  CHROMA 

.36  10 

.69  100 

1.07  0  2.71  ± 

1.42  0 


100 

100 

.00  0  1.41  ± 

60 


100 

100 

.53  100*088.70  ±*0.32 

100 


MINUS  CHROMA 

.36 

SO 

100 

.69 

100 

100 

1.07 

10  .43  ± 

.16  0 

1.42 

0 

50 

PLUS  HUE 

.18 

90 

100 

.38 

100 

100 

.57 

100  .72  ± 

.31  100 

.75 

0 

0 

MINUS  HUE 

.18 

100 

100 

.38 

100 

100 

.57 

100  .67  ± 

.04  100 

.75 

10 

10 

PLUS  LIGHT 

.50 

100 

100 

1.00 

10 

100 

1.50 

10  .80  ± 

.15  0 

2.01 

0 

0 

MINUS  LIGHT 

.50 

100 

90 

1.00 

10 

90 

1.50 

0  .77  ±  .13 

0 

9.01 

0 

0 

100 

100 

1.24  ±  .32  100*088.70  ±*0.32 

100 


100 

100 

.65  ±  .04  100*472.55  ±*2.88 

100 


.67  ±  .04 


90 

100 

100 

80 


21*81  ±*0.04 


1.22  t  .09 


100 

100 

100 

10 


1.77  ±  .12 


100 

100 

1.15  ±  .23  100 

0 


1.74  ±  .12 
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Table  7 


Raw  Results  for  Observer  G  on  Dark  Blue,  Together  with 
Calculated  50%  Acceptance  Limits 


PCT.  BOUNDARY 
PASS  4  STD.  DEV. 


80 

100 

10  .91  4  .24 


40 


100 

90 

0 

100 


1.82  4  1.45 


100 

100 

90*443.87  4*5 . 12 


100 


100 

90 

100 

30 


.68  1  .16 


90 

90 

20 

10 


1.15  4  .1* 


100 

100 

60 


1.55  4  .14 


Table  8.  Supmary  of  Tolerance  Results 


Olive  Green 

+  Chroma 

-  Chroma 

+  Hue 

-  Hue 

+  lightness 

-  lightness 

A 

0.78 

0.83 

1.76 

B 

1.05 

0.46 

0.53 

1.64 

1.78 

c 

1.47 

0.71 

0.53 

2.29 

D 

1.38 

0.92 

0.49 

0.55 

1.46 

1.58 

E 

1.13 

0.85 

0.55 

0.55 

1.50 

1.42 

F 

1.42 

1.21 

0.83 

0.60 

+  St. 

D. 

1.13  t 

0.28 

0.60  ± 

0.13 

1.68  £  0. 

28 

Tan 

•>  Chroma 

-  Chroma 

+  Hue 

-  Hue 

+  lightness 

-  Lightness 

A 

0.8S 

0.73 

0.25 

1.41 

B 

0.85 

0.77 

0.43 

0.29 

1.27 

1.50 

C 

0.98 

0.80 

0.30 

1.09 

1.15 

D 

1.04 

1.22 

0.54 

0.59 

1.07 

0.93 

E 

1.92 

0.57 

0.55 

2.55 

2.64 

+  St. 

D. 

1.03  t 

0.37 

0.44  t 

0.14 

1.51  ±  0 

.64 

Dark  Blue 

+  Chroma 

-  Chroma 

♦  Hue 

-  Hue 

+  lightness 

-  Lightness 

A 

.59 

1.04 

B 

.61 

1.29 

1.17 

G 

.93 

.56 

.62 

D 

.67 

.80 

.77 

E 

1.24 

.65 

.67 

1.22 

1.15 

G 

.91 

.68 

1.15 

1.55 

■+•  St- 

n. 

0.92  ± 

0.27 

0.66  t  0.03 

1.03  t  0 

.31 

1.  Differences  between  plus  and_  minus  directions..  There  is  no  8  ^ 

difference  between  the  plus  hue  and  minus  hue  tole^es  or  between  the  p  us 
lightness  and  minus  lightness  tolerances.  The  difference  between  the^lus^ 
chroma  and  minus  chroma  tolerances,  however,  is  debat  •  “  di£ference 

Student's  t  test  to  the  unpaired  data,  there  is  ™  and 

between  plus  and  minus  chroma.  If,  however,  we  pair  t  P  larger 

minus  chroma  data  by  observer,  we  see  that  three  observers  h^aurge^ 
plus  chroma  tolerance  than  a  minus  chroma  tolerance  for  “t,sSt  test 

whereas  no  observer  had  the  reverse.  Application  of  the 

to  the  paired  data  showed  significance  for  this  color,  u  observers 

even  the  paired  data  did  not  show  a  significant  difference;  three  observers 

favored  plus  chroma  and  one  favored  minus  chroma.  e  ec ections  We 
that  no  difference  existed  between  any  plus  and  minus  c.icuiattng 

therefore  avar.gud  all  the  values,  plus  and  *i£,1£1^“ln  T.JU  8. 

averages  and  standard  deviations,  and  these  latter 

2 .  Differences  Between  Hue,.  Ch  r ratios 

1.  strong"" confirmation  of  the  ’  o  U«  g  « “  trifle, 

in  the  oases  of  the  olive  green  and  the  tan  For  the  o U«  8  ^ 

«;s"ndirr:;.cir“t?  »eu:.:"h:t8it  ls  ^rr  „  .-.mu-. •»* 

Ve  must  work  out  individual  tolerances  for  each  color. 

The  absolute  value,  for  the  hue  ClK"' 

hood  of  half  a  CIELAB  unit,  t  ose  or  c  r  one_half  CIELAB  units  for 

unit.  The  lightness  tolerances  are  ne«  JjJ  °nit  for  the  dark  blue, 

the  olive  green  and  the  tan,  and  near  one  CIELAB  unit  tor 

,  Differences  between  Observers  from  Industrj^^ndjrgm^^ 

Service.  We  might  possibly  expect  that  the  observers  from  n  ustry  wou 

^r^trca^f™  rai:rrh^ether:in;Tfhtere8iS  no 
sign! f leant *" difference  between  the  results  of  the 

i^b^Von  fhrtrnand°theeTarrk9ble,  but  his  results  for  these  colors 
were  not  included  in  Table  8. 

IV.  ACCEPTABILITY  EQUATIONS 

ellipsoid  determines  whether  a  given  sample  does  or  does  not  pass. 


to 


This  equation  takes  the  following  form 

2 


A  A  = 


fu  (<la*)  + 


2g 


12 


Aa  Ab  +  g22 


(Ab*)2  +g33  (AL") 


(1) 
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where  A  A  is  an  acceptability  figure,  so  scaled  that  normally  if  A  A  is  less 
than  1,  the  sample  is  accepted;  ifA  A  is  greater  than  1,  it  is  rejected  (see 
further  discussion  on  this  point  below).  The  quantities  A  a*,  A b*,  A  L*  are 
sample  minus  standard  in  CIELAB  coordinates. 


’22 


Appendix  G  explains  the  derivation  and  meaning  of  the  acceptability 
equation.  The  coefficients  g^,  2g.2»  8-i 
following  equations  where  a*^  and  b*Q 
values  of  the  standard,  c  is  the  chroma  tolerance,  h  is  the  hue  tolerance, 
and  v  (for  value)  is  the  lightness  tolerance: 


and  g~,  are  given  by  the 
sent  tne  CIELAB  a*  and  b* 


0 

=  tan  1 

(>; 

(2) 

*11 

(  2 
=  I  cos 

e/c2) 

1  +  ^sin2  0/h2^  , 

(3) 

2*12 

\ 

=  2  sin 

0  cos 

0  (l/c2)  -  (l/h2)j  , 

(4) 

*22 

=  (sin2 

0/c2 

^  +  ^cos2  0/h2^  , 

(5) 

*33 

=  1/v2. 

(6) 

To  illustrate  how  these  equations  are  used,  let  us  apply  them  to  the 
olive  green.  The  olive  green  standard  has  the  CIELAB  coordinates  L*n  = 
31.71,  a*  =  3.76,b*  =  9.31  (see  Table  9).  We  have  0  =  tan  ”  (9.31/- 

3.76)  =  68.01°.  U 

This  gives  g^  =  [cos2  (-68.01) /l.  132]  +  [sin2  (-68.01) /0. 602]  =2.50. 
Similarly.  2gJ2  =  1.39;  g22  =  1.06,  g33  =  0.354.  We  therefore  have 


1.06 


(a*  +  3.76)  2  +  1.39  ^a*  +  3.76)  (b*  - 
(b*  -  9.3l)  2  +  0.354  (L*  -  31. 7l)  2  , 


9.31  + 


(7) 


where  a*,  b*  and  L*  are  the  CIELAB  coordinates  of  the  sample  being  judged. 
This  equation  will  have  to  be  modified  each  time  the  standard  is  changed. 
The  values  for  c,  h  and  v,  of  course,  were  taken  from  Table  8  and  would 
remain  constant  irrespective  of  the  standard. 
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The  equations  for  the  other  colors  are  derived  in  the  same  way.  The 
CIELAB  coordinates  for  the  standards  are  given  in  Table  9.  We  have  data 
for  two  different  sets  of  limit  samples  that  were  used  with  the  tan;  each 
set  had  its  own  standard.  Following  are  the  equations  to  be  used  with 
each  of  the  tan  standards  and  the  equation  for  the  dark  blue. 


For  the  standard  of  tan  set  no.  Is 
A  A  = 


5.09  (a*  -  2.05)  2  -  1. 09 ^a*  -  2.05^  (b~  -  15.64) 


1.01  (b*  -  15.64)  2  +  0.439  ^L*  -  57.48) 

For  the  standard  of  tan  set  no.  2: 

A  A  =  |  5.06  (a*  -  2.47)  2  -  1.35  (e*  -  2.47)  (b*  -  15.09)  + 
1.05  (b*-  15.09)  2  +  0.439  ^L*-  56.43)  2  , 

For  the  dark  blue: 

AA  =  !  2.30  (a*  -  0.205)  2  +  0.041  (a*  -  0.205)  (b*  +  11. 12) 
1.18  (b*  +  11.12)  2  +  0.943  (l*  -  22.98)  2  , 


(8) 


(9) 


(10) 


We  tried  these  equations  on  the  limit  samples  for  each  color.  Results 
are  shown  in  Table  9.  We  see  that  all  of  the  limit  samples  for  the  olive 
green  and  most  of  those  for  tan  set  no.  1  would  fail  the  acceptability  tes 
as  given.  The  average  A  A  values  are  1.65  for  the  olive  green  and  1.63  for 
tan  set  no.  1.  For  tan  set  no.  2  th-  average  A  A  drops  to  1.21,  whereas 
for  the  dark  blue  the  average  A  A  is  0.76,  with  ail  but  one  of  the  limit 
samples  passing  the  test. 


It  is  widely  believed  in  industry  that  for  commercial  customers  the 
range  of  limit  samples  is  influenced  by  the  state  of  the  marketplace  at 
the  time  of  choice.  If  goods  are  needed  urgently,  the  limit  sample  mighi 
be  chosen  so  as  to  reflect  wider  tolerances;  if  there  is  enough  material 
and  the  concern  is  more  with  high  quality,  the  reverse  might  be  true. 
This  is  rarely  true  of  military  procurement;  other  factors  are  far  more 
important.  Industrial  experience  in  matching  a  standard  may  result  in 
tightening  the  range,  as  may  be  seen  by  comparing  tan  set  no.  1  with 
tan  set  no.  2;  the  tolerance  is  considerably  tighter  in  the  latter  set. 
It  is  tighter  still  in  the  dark  blue  set  because  the  end  use  of  the 
fabric  is  mote  critical. 
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Table  «  Application  of  Tolerance  Equations  to  limit  Samples 


DESIG. 


* 

L 


* 


a 


* 

b 


A  A 


Olive  Green 

STD 

FY 

FS 

FG 

FB 

FR 

TG 

TY 

TB 

TS 


31.71 
32.29 

31.72 
31.87 
32.35 
32.10 
32.23 
32.43 
33.79 
32.65 


-3.76 

-5.29 

-4.92 

-3.85 

-4.28 

-4.65 

-3.74 

-4.26 

-4.87 

-4.89 


9.31 

9.76 

9.45 

8.03 

8.51 

9.57 

8.07 

8.78 

8.80 

9.75 

Av.  ±  St.  D. 


2.29 
1.78 
1.39 
1.44 
1.34 

1.30 
1.21 
2.38 
1.74 

1.65  ± 


Tan.  Set  No.  1 

STD 

FS 

FY 

FG 

FR 

FB 

TS 

TY 

TR 

TB 


57.48 

56.47 

57.26 

57.54 

58.52 

58.00 

58.96 

58.50 

58.42 

57.94 


2.05 

3.31 

2.92 

2.49 

3.41 

2.35 

2.21 

2.53 

2.67 

2.14 


15.64 

16.38 

16.91 
16.36 

15.91 
15.27 
15.34 
16.74 
16.30 
15.27 

Av.  ±  St.  D. 


2.84 

2.07 

1.08 

3.09 

0.91 

1.11 

1.51 

1.53 

0.56 

1.63 


Tan.  Set.  No.  2 

STD 

FS 

FR 

TY 

FY 

FG 

TG 

TB 

FB 

TS 

TR 


56.43 

56.06 

56.89 

56.92 

55.66 
55.97 
56.34 
57.23 
56.84 
57.83 

57.67 


2.47 

2.80 

3.21 

1.86 

2.24 

1.82 

1.93 

2.57 

2.82 

2.56 

2.88 


15.09 

15.07 

15.23 

13.60 

15.60 
15.12 
15.21 

14.71 
14.50 
14.63 

15.72 


Av.  ±  St.  D. 


0.79 

1.66 

1.64 

0.98 

1.50 

1.26 

0.73 

1.16 

1.09 

1.26 

1.21 
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0.43 


0.87 


±  C. 33 


DES1G. 

Dark  Blue 

STD  1  * 

STD  2'  * 

FG  1 

FG  2 

TG 

TS 

TB 

FB 

FR 

FS 

TR 


Table  9 

(cont'd 

* 

* 

L 

a 

22.90 

0.40 

23.06 

0.01 

20.93 

0.16 

22.18 

-0.06 

23.09 

0.05 

22.53 

0.23 

23.56 

-0.03 

22.70 

-0.25 

22.34 

-0.16 

22.59 

0.04 

22.85 

0.13 

Av. 


* 

b 

A  A 

-11.49 

-10.75 

-10.81 

2.02 

-10.79 

0.94 

-11.48 

0.47 

-10.89 

0.51 

-11.16 

0.67 

-11.25 

0.76 

-10.92 

0.86 

-10.96 

0.48 

-11.12 

0.17 

±  St.  D. 

0.76  ±  0.53 

*  Two  readings  on  same 


standard;  values  were  averaged 


The  tolerances  for  the  calculation  method  presented  here  can  be  varied 
by  merely  changing  the  limit  for  <dA  at  which  we  are  willing  to  accept  samples 
upward  or  downward  from  1.0.  For  example,  if  the  A  A  limit  is  set  at  1.6 
instead  of  1.0*  we  would  have  approximately  the  tolerances  that  obtain  for 
the  olive  green  limit  samples.  If,  on  the  other  hand,  we  set  the  A  A  limit 
at  0.8,  we  would  have  the  more  stringent  tolerances  that  hold  for  the  dark 
blue  limit  samples.  A  tolerance  limit  of  1.0  represents  average  practice, 
as  determined  by  the  judgments  of  unbiased  inspectors  obtained  as  described 
in  the  sections  above.  We  should  remember  that  any  change  in  tolerance  limits 
as--just  described  will  still  maintain  the  relative  importance  of  the  hue, 
chroma  and  lightness  directions;  we  are  merely  changing  the  overall  tolerance. 

V.  CONCLUSIONS  AND  RE 00MMEN DATION S 

The  methods  described  in  this  report  for  selecting  samples  for  inspection, 
for  presenting  the  samples  to  the  observers,  and  for  scoring  the  results  and 
setting  up  the  acceptability  equations  appear  to  be  reliable.  For  the  present, 
we  recommend  that  every  color  for  which  tolerances  are  desired  by  the  Army  be 
handled  In  the  same  way  as  we  have  illustrated  for  ".he  olive  green,  the  tan 
and  the  dark  blue.  The  procedure  allows  for  flexibility  in  setting  up  the 
tolerances  by  simply  changing  the  A  A  acceptability  limit  upward  or  down¬ 
ward  from  1.0.  This  choice  can  be  made  anew  each  time  any  change  in  pro¬ 
curement  requirements  occurs. 

It  may  eventually  happen  that  enough  history  and  experience  will 
accumulate  to  make  it  unnecessary  to  obtain  observers*  judgments  in  setting 
up  equations  for  each  new  color.  We  may  see  enough  regularity  in  the 
ellipsoids  to  make  it  unnecessary  to  determine  rr.ai:y  more.  Considerably 
more  data  like  those  in  Table  8  would  be  needed,  however,  before  we  could 
make  any  such  generalizations. 
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Appendix  A.  Computer  pi’ogram  for  study  of  relative  positions  of 
standard  and  limit  samples  (Program  RAMSLEY) 


This  program  constructs  all  possible  tetrahedra  from  the  nine  points  in 


color  space  corresponding  to  the  standard  and  the  eight  limit  samples  (seven  limit 

''q 

samples  in  one  case).  There  are  ' 


=  126  such  tetrahedra  that  can  be  constructed 
v*/ 

from  the  nine  points. 

The  program  takes  each  of  the  nine  points  in  turn,  and  determines  whether 


it  falls  inside  any  one  of  the 


/8\  . 
W 


70  tetrahedra  that  are  not  constructed  from  the 
point  in  question.  If  the  point  falls  inside  at  least  one  of  the  tetrahedra,  it  is 
obviously  on  the  inside  of  the  geometric  figure  determined  by  the  nine  points;  if 
not,  it  is  on  the  boundary.  (We  define  the  "geometric  figure  determined  by  the 
nine  points"  as  the  union  of  all  126  tetrahedra.) 

The  determination  of  whether  a  point  falls  inside  a  tetrahedron  is  made  as 
follows:  Call  the  four  corners  of  the  tetrahedron  A,  B,  C  and  D;  the  point  in 
question  is  called  P.  Choose  arbitrarily  the  D  corner  of  the  tetrahedron  and 
construct  vectors  to  the  other  three  corners  (see  Figure  5).  Also  construct 
a  vector  from  point  D  to  the  point  in  question  (dotted  line).  Call  the  vectors  to 
the  three  corners  DA,  DB  and  DC;  the  vector  to  the  point  in  question  is  DP. 

Now  the  point  is  inside  the  tetrahedron  if 

_ i  _ \  i  ^ 

pDA  +  qDB  +  rDC  =  DP;  p  >  0; 

q  >  0; 
r  >  0; 

1-p-q-r  >  0. 


Expressing  the  vector  equation  in  algebraic  notation,  if  L A  indicates  the  L 
value  of  point  A,  etc. ,  we  have  the  following  matrix  equation: 
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The  computer  program  forms  all  126  matrices  shown  on  the  left  above, 
and  inverts  them.  It  then  chooses  the  70  inverted  matrices  to  be  used  for  each 
point,  multiplies  them  in  turn  into  the  appropriate  right-hand  side  vector,  and 
calculates  p,  q  and  r.  If  the  inequalities  shown  above  for  p,  q  and  r  are  satisfied 
in  at  least  one  of  the  70  cases ,  the  point  in  question  is  inside  the  geometric  figure 
formed  by  the  nine  points. 


The  computer  program  follows. 


noo  o o  o  o  ooo  o  oo  o  o  o  o  o  n  o 


PROGRAM  RAMSLEY  (INPUT,  OUTPUT •  TAPE1 *  TAPE2) 

C 

C  FOR  A  SERIES  OF  12  COLORS*  READS  THE 

SPECTROPHQTOHETRIG  CURVES  OF  THE  STANDARDS 
AND  THE  LIMIT  STANDARDS*  ANO  FOR  EACH  OF 
THESE  DETERMINES  WHETHER  ITS  POINT  IN 
CIELA8  COLOR  SPACE  LIES  ON  THE  BOUNDARY  OF 
THE  THREE-DIMENSIONAL  FIGURE  DETERMINED 
8Y  ALL  THE  POINTS. 

DIMENSION  NAME(6*  12)*  LIMNUMI9,  12),  A(9),  8(9), 

♦I0(*,  126),  RHS  ( 3) ,  C(*l  «  LA8ELC9)  •  WORK(S),  LIM(2,  9) 

REAL  L  ( 9)  *  MAT  (3*  3),  MIN  V  ( 3*  3),  MINVSI3,  3*  126) 

OATA  NANE/8H0LIVE  DR*  8HA3  7,  UN, 8HTREATED  .SHCOTTON  O* 
♦SHUCK,  d. 2  »8H5  OZ.  «  dHT  AN  46.  *  8HCQTT  ON  P.BHOPLIN,  4, 

♦  6H  OZ.  ,8H  ,8H  * 8HT  AN  M-1,,8H  CL.  POL, 

♦8HY/W00L  T  *  8HR0P.  9  Q,8HZ.  ,8H  ,8HARMY  aRt* 

♦  8HEN  4*. ,  W,8HL  GAB.*  *8H130Z.  *  8H  *8H  * 

♦oHOlIVE  GR.8HEEN  1J7  *  *  dH  NYLON/C, BHTN  (50/5, 8H0)  P3PLI, 

♦  8HN,  5  OZ .  *  8HBLUE  150.8HTROP.  WL,8H.  10  0Z.,8H  * 

♦  8H  ,8H  , 3HARM Y  GRE, 8HEN  3*4,  » 8HP0LY ESTE* 

♦  8HR/N00L  G ,8H A3 AR DIN E * BH  9  OZ.  *  8HARMY  GRE,8HEN  344,  » 

♦dHPOLYESTE , 8HR/W00l  T,  8HR0PICAL, * 8H  10  OZ .  » 8 HA RM Y  GRE, 

♦8HEN  44*  W, 8HL.  SERGE, 8H  15  OZ.  , 8H  ,8H  , 

♦  8H0LIVE  GR, SHEEN  108,, dH  WOOL/NY,  8HL0N*  16  ,8H0Z.  SHIR* 

♦8HTING  , 8HOLIVE  GR, SHEEN  107, ,SH  CTN.  SA.SHTEEN,  d., 

♦  8H2  OZ.  ,8H  ,  SH3LUE  151.8H*  TROP.  ,8HML.,  10  , 

♦  8HOZ.  * 8H  *  SH  / 

DATA  LIMNUM/1,2,3,4,  5,7 ,8*9,10,1,2,5,6,7, 8, 10*11, 3, 

♦  It  2,  3, 4, 5, 7*8, 9*10*1*2*6, 4*5*7,  11  *9, 10,  1,  2*  3,  *«  5, 7  *8,9,10, 

♦  1,2, 6, 4, 5*7, 11, 9*10,  1*2, 6*4*5, 7, 11*9, 10, 1*2*6, 4, 5, 7  *11,9,10, 

♦  1,6,4*5,7,2,11,9,10,1,2,3,4*5*7,8*9,10,1,2*3,5,  *♦, 7, 8*9, 10, 

♦  1,  2,6,4,5,7,11,9,10/ 

START  LARGE  00  LOOP  BASED  ON  12  SAMPLES 
REWIND  1 

00  500  NS TO  =  1,  12 

START  00  LOOP  8ASE0  ON  NINE  STANDARDS  ANO  LIMIT 
SAHPLES  (EIGHT  IN  THE  CASE  OF  SAMPLE  NO.  2) 

NP  *  9 

IF  (NSTO  *  EQ  .  2)  NP  =  8 
00  190  M  -  l*  HP 
READ  (1*  1QG1)  L(M),  AIM),  B(M) 

ASSIGN  NAMES  TO  STANOAROS  AND  LIMIT  STANDARDS. 

INDEX  *  LIHNUH4H,  NSTO) 

GO  TO  (171,  172,  173,  174,  175,  176,  177,  178,  179,  180, 
♦181),  INDEX 

l?t  Lir*(i,  M)  =  ttHSTANOARQ 

LIM(2,  M)  -  SH 

GO  TO  199 

172  LINI1,  HI  «  8HTHIN  STA 
LI M ( 2 ,  HI  -  8HNOARO 
GO  TO  190 

173  LIHll,  M)  =  8HTHIN  Y EL 

LIMI2,  Ml  i  8HL0W 


0010  00 
001010 
0  Q 10  20 
0010  30 
001046 
0010  50 
001060 
0010  70 
001080 
0  010  90 
001100 
001110 
001140 
001150 
001160 
001170 
001180 
001190 
001200 
001210 
001220 
001230 
0012*0 
001250 
001260 
001270 
001280 
001290 
001300 
001310 
001320 
001330 
0013*0 
001350 
001360 
001370 
001389 
001490 
001500 
001505 
001510 
00152Q 
0015  30 
001540 
Q0155A 
001560 
001570 
001655 
001660 
001910 
001920 
001930 
001940 
0  019  50 
001960 
001970 
SO  19 80 
001990 
0020  OS 
0020 10 
0020  20 
0020  30 
0020*0 
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GO  TO  190 

174  LIHil,  Ml  *  dHTHIN  GRP 

LIM( 2,  Ml  =  8HEN 

GO  TO  190 

175  LIHil,  HI  =  8HTHIN  RED 

LIM( 2,  HI  =  dH 

GO  TO  190 

176  LIHil,  HI  =  8HTHIN  8LJ 
Li  Hi  2 »  HI  =  8HE 

GO  TO  190 

177  LIHil,  HI  =  8HFULL  5 TA 

LIHI2,  HI  =  8HN0AR0 

GO  TO  190 

17d  LIHil*  HI  =  dHFULL  YEw 

L IHi 2*  HI  s  8HL0H 

GO  TO  190 

179  LIHil,  Ml  =  8HFULL  GRE 

LlMi  2*  HI  =  dHEN 

GO  TO  190 

180  LIHil,  Ml  =  8HFULL  RED 
Li  Hi  2,  Ml  =  8H 

GO  TO  190 

181  LIHil,  HI  =  8HFULL  3LO 

LlHi 2,  HI  =  8 HE 

190  CONTINUE 

CREATE  MATRICES  REPRESENTING  TETRAHEDRA  FOR  ALL 
THE  POSSIBLE  COMBINATIONS  OF  THE  POINTS  TAKEN 
FOUR  AT  A  TIME 

NHAT  =  0 

NPM1  =  NP  -  1 

NPM2  s  NP  -  2 

NPM3  *  NP  -  3 

DO  200  J1  =  1,  NPH3 

JlP  =  Jl  ♦  1 

DO  200  J2  =  JlP,  NPH2 

J2P  *  J2  *  i 

DO  200  J3  =  J2P,  NPH1 

j3P  :  J3  U 

OQ  200  J4  =  J3P,  NP 

NHAT  *  NHAT  ♦  1 

10(1,  NHAT I  s  ji 

IOi2,  NHAT I  =  J2 

10(3,  NHAT I  =  J3 

10(4,  NHAT  I  =  J* 


HATH, 

11 

S 

LCJ1I 

m 

LC  J4I 

MAT  Cl, 

21 

= 

L  <J2 1 

- 

L  (  J4 1 

HAT ( 1 , 

31 

s: 

L(J3I 

«* 

Li  J4I 

MAT (2, 

11 

at 

A  ( Jll 

• 

A  (  J4  I 

MAT ( 2, 

21 

= 

A  (J2I 

- 

AC  J4I 

MAT ( 2, 

31 

3 

A  (j3  1 

m 

A{  J4I 

MAT (3, 

11 

S 

3  (Jll 

• 

S(  J4I 

MAT (3, 

21 

3 

S(J2I 

«• 

B(  J4I 

MAT43, 

31 

= 

8  ( J3I 

• 

B(  J4I 

INVERT  MATRICES  AND  STORE  THEM  IN  MEMORY. 

CALL  UNVlF (MAT ,  3,  3,  KINV,  0,  HORK,  IERI 
OO  200  J  =  1,  3 
OO  200  K  =  1,  3 

200  HI  NV  S  ( J,  K,  NHAT I  *  HXNViJ,  Kl 

TRY  EACH  POINT  SEPARATELY  KITH  EACH  OF  THE 
INVERTED  MATRICES,  TO  SEE  IF  IT  FALLS  HlTMlN 
THE  TETRAHEDRON  DEFINED  BY  THF  FOUR  POINTS  THAT 


00  24  50 
002060 
0  0  20  70 
Q  02080 
002090 
002100 
002110 
002120 
00  21 30 
002135 
0  02140 
002150 
00216Q 
002170 
002180 
002190 
CO  22  00 
002210 
002220 
002230 
0022*0 
002250 
002260 
002270 
002310 
002320 
002330 
002340 
002350 
002360 
002370 
002380 
002390 
0024  00 
002410 
00  2420 
002430 
002440 
0  0  2450 
002460 
002470 
0Q2480 
00  2490 
00  25  0C 
002510 
002520 
002530 
002540 
002550 
002560 
002570 
002580 
002590 
002600 
002610 
002620 
0026 30 
002640 
002650 
0  02660 
002670 
002680 
C0269C 
002700 
002710 
c  fl?7  ?n 


o  o  o 


C 

C 

C 


WENT  INTO  THE  MATRIX--  l(  SO*  LABEL  THE  POINT 
AS  BEING  IN.ilOE  THE  BOUNDARY. 


00  300  J  *  It  9 
DO  270  NHAT  *  i»  12& 

IF  (J  .  EQ  .  IDClt  NNAT)  . 

♦  J  •  EQ  •  10(3*  NMATI  .  OR 
JJ  =  10(4*  NHAT ) 

RHSl 11  s  LIJ)  -  LUUI 
RHS(2)  =  A(J)  -  A  (JJ ) 
RHS(3»  =  3  (J )  -  B(JJ» 

00  250  K  =  1*  3 


OR  .  J  • 

.  J  .  £3 


£Q  .  10(2*  NHAT)  .  OR  • 

.  10(4*  NHAT) )  GO  TO  270 


coo  =  o. 

DO  250  KK  = 
250  C ( K)  =  COO 
C<4)  =  1.  - 
IF  (C(ll  • 

♦  C (3)  .  GT 


1,  3 

♦  HINVS(K*  KK*  NHAT) 


*  RHS(KK) 


C(i)  -  C (2)  -  C (3) 

GT  .  0.  AN3  .  C<2)  . 
.  0.  .  AND  .  C<4)  .  GT 


QT  •  0.  •  AND  • 

.  0.)  GO  TO  280 


270  CONTINUE 

LA3ELCJ)  *  2 
GO  TO  300 


280  LABEL (J)  =  1 
300  CONTINUE 


PRINT  RESULTS 

PRINT  1002*  (NAHE(J«  M$TO)«  J  =  1«  6) 

PRINT  1008 
00  430  J  =  1*  NP 
INOEX  *  LABEL (J) 

GO  TO  (310.  320).  INOEX 

310  PRINT  1003,  LIH(1.  J).  LIH(2.  J>»  LCJ).  AUl,  BCJ 

320  PRINT  1004,  LIH  ( 1,  J)»  LINI2,  J) .  L(J),  ACJ).  B(J) 

<*00  CONTINUE  „„„  _n  *  nB 

IF  ( NSTQ  .  EQ  .  4  .  OR  .  NSTO  .  EQ  .  8  •  OR  . 

f  NSTD  .  EQ  .  12)  PRINT  1810 
IF  (NSTD  .  EQ  .  2)  PRINT  1006 

500  CONTINUE 
STOP 


FORMAT  STATEHENTS 
1001  FORMAT t3Fd. 2) 

1302  FORMAT C///6X*  6A8)  BOUNDARY*) 

1003  FORMAT (IX,  2A8,  F7.2,  2F8.2*  3X,  *NOT  ON  BOUNDARY 

ilti  FORHATUX,  2A«.  F7.2.  2FA.2,  JX,  »ON  aOUNO*R»*> 

1006  FORMAT C 1H  )  M1  _  -A-  #a*. 

1008  FORMAT C/3X,  * SAMPLE* «  12X,  L  *  7X,  A  ,  7  . 

♦  7X.  *POSI TION*/ ) 

1010  FORMAT (/) 

END 


0  027  30 
002740 
0027  50 
0  0  27  60 
002770 
00  27  80 
002790 
0  0  28  00 
002810 
00  28  20 
0028  30 
0028  40 
002850 
0  0  2860 
002870 
0  028  80 
0028  90 
0  0  29  00 
002910 
0  0  29  20 
0029  30 
002940 
002950 
002960 
002970 
002980 
002990 
003000 
0030 10 
00  30  20 
0030  30 
003040 
00  30  50 
003060 
003105 
003110 
003115 
003116 
003120 
00  3130 
0  0  31**Q 
003150 
003160 
003170 
003180 
003190 
003200 
003220 
0Q32«*0 
003250 
003270 
003280 
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Appendix  B.  Computer  program  for  determination  of  minimum  non-tilting 
ellipsoid  containing  standard  and  limit  samples  (Programs 
ELPSOID  and  ELPSOD2) 


A  non-tilting  ellipsoid  (in  the  sense  described  ir.  the  main  scclion  of  this 
report)  is  completely  determined  by  seven  parameters:  the  three  position  coordinates 
of  the  center,  the  three  semiaxis  lengths,  and  the  angle  between  the  major  horizontal 

j|c 

axis  and  the  a  axis.  For  Program  ELPSOID,  all  seven  parameters  are  adjustable. 


For  Program  ELPSOD2,  only  six  parameters  are  adjustable;  the  axis  angle  is 

constrained  by  the  fact  that  the  major  horizontal  axis  lies  along  a  line  connecting 

* 

the  center  of  the  ellipsoid  with  the  L  axis. 

Consider  a  non-tilting  ellipsoid  located  somewhere  in  the  a  -b  -L  coordinate 
system.  Set  up  another  coordinate  system  to  coincide  with  the  three  axes  of  the 
ellipsoid.  The  equation  of  the  ellipsoid  in  this  second  coordinate  system  is 


2  2  2 

+  5_ 

2  2  2 

p  q  r 


(B-l) 


where  x ,  y  and  z  are  the  coordinates  of  any  point  on  the  ellipsoid  relative  to  the 
second  coordinate  system,  and  p,  q  and  r  are  the  semiaxis  lengths.  We  consider 
that  x  and  y  are  the  axes  parallel  to  the  a  -b  p^ane  and  z  is  the  axis  that  is 

Jjf 

parallel  to  the  L  axis.  The  major  horizontal  semiaxis,  of  length  2p ,  lies  along  the 
x  axis. 

* 

If  9  is  the  angle  going  from  the  positive  branch  of  the  a  axis  counter- 

*  *  * 

clockwise  to  the  positive  branch  of  the  x  axis,  and  if  ac,  b£  and  Lc  are  the 

£  3|c  j|c 

coordinates  of  the  center  of  the  ellipsoid  in  the  a  -b  -L  system,  then  the  following 
transformation  equations  apply: 

4c  4*  4*  jj* 

x  =  (a  -  a  )  cos  9  +  (b  -  b  )  sin  e, 
c  c 


]|c  4  jf;  ]j( 

y  =  -(a  -  a  )  sin  9  +  (b  -  b  )  cos  0, 
c  c 


(B-2) 


z  =  L-L 
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Substituting  these  equations  into  the  equation  of  the  ellipsoid  given  above  and 
collecting  terms,  we  have 


This  is  the  equation  for  the  ellipsoid  in  the  a  -b  -L  coordinate  system 

incorporating  the  seven  parameters  mentioned  above.  The  position  parameters  for 

the  center  are  a  ,  b  and  L  ;  the  semiaxis  lengths  are  p,  q  and  r;  the  angle 
c  c  c 

* 

between  the  major  horizontal  axis  and  the  a  axis  is  0. 

Program  ELPSOID  works  in  the  following  way:  The  se  n  parameters  are 
systematically  varied,  one  at  a  time.  After  each  change  in  the  parameter.-,  they 
are  insei-ted  in  Equation  B-3  to  give  what  may  be  called  a  reference  ellipsoid.  All 
the  points  (a  ,  b*.  L  )  that  satisfy  Equation  B-3  fall  on  the  surface  of  this 
reference  ellipsoid. 

For  each  parameter  change  and  associated  reference  ellipsoid,  the  program 

then  enters  a  subroutine  called  SIZEIt.  This  subroutine  goes  through  the  nine 

points  for  the  standard  and  the  limit  samples,  and  substitutes  the  coordinates  of 

each  point  at  a  time  into  the  equation  for  the  reference  ellipsoid.  Consider  one  of 

these  points.  Since  this  point  is  probably  not  on  the  reference  ellipsoid,  the  right 

2 

hand  side  of  the  equation  will  not  be  unity  but  some  value  that  we  will  call  S  . 

2  t  i  * 

Dividing  Equation  B-3  through  by  S  ,  and  setting  p  =  pS ,  q  =  qS  and  r  =  rS, 

t  t 

we  obtain  a  new  equation  just  like  Equation  B-3  but  with  p  replacing  p,  q 

I 

replacing  q,  and  r  replacing  r,  and  with  unity  back  on  the  right  hand  side. 
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This  new  equation  defines  a  new  ellipsoid  that  is  concentric  with  the  reference 
ellipsoid ,  that  has  the  point  in  question  on  its  surface  and  that  has  the  same 
ratios  of  axis  lengths  as  those  of  the  reference  ellipsoid,  since  the  new  axes  are 
merely  the  old  axes  multiplied  by  S.  The  volume  of  an  ellipsoid  is  equal  to  4/3  ir 
multiplied  by  the  product  of  the  three  semiaxes;  the  new  ellipsoid  therefore  has  a 

3 

volume  equal  to  4/3  tt  times  the  product  of  the  p,  q  and  r  parameters  times  S  . 

The  subroutine  goes  through  all  .nine  points  in  this  way,  and  computes 
nine  corresponding  ellipsoids ,  all  concentric  with  and  similar  (in  the  geometric 
sense)  to  the  reference  ellipsoid.  It  then  returns  to  the  main  program  the  volume 
of  the  largest  of  these  concentric  ellipsoids ,  since  that  is  the  one  that  contains 
all  the  points.  If  this  volume  is  smaller  than  the  smallest  volume  previously 
determined  for  other  parameter  variations,  the  current  variation  wa3  favorable; 
otherwise  not.  The  systematic  parameter  variation  continues  until  no  parameter 
change  can  produce  a  smaller  ellipsoid  returned  by  the  subroutine.  The  program 
then  normalizes  the  current  parameters  by  multiplying  the  p ,  q  and  r  values  by 
the  S  value  for  the  ellipsoid  that  was  chosen,  and  finally  prints  out  the  ellipsoid 
characteristics. 

In  Program  ELPSOD2,  the  angle  0  is  always  fixed  by 


g33 
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PROGRAM  ELPS0I0  (INPUT,  OUTPUT,  TAPE3) 

FOR  A  SERIES  OF  POINTS  IN  CIELA3  COLOR  SPACE, 

DETERMINES  THE  SMALLEST  NON-TILTING  ELLIPSOID 
THAT  CONTAINS  THE  POINTS.  THE  ELLIPSOID  MAY 
BE  INCLINED  WITH  RESPECT  TO  THE  HORIZONTAL  AXIS, 

BUT  MAY  NOT  TILT  EITHER  FORWARD  OR  SIOEWAYS. 

DIMENSION  PA  R  (7  )  ,  H0WFARI9),  OELTA(7) ,  014(4),  NAME (6 ,  12), 
♦SIGN (2) 

COMMON  CIE (3,  9),  NP 
LOGICAL  CHANGE,  ACTION 

DATA  DELTA/C. 1,  0.1,  0.1,  0.1,  0.1,  0.1,  0. 017453293/, 
♦014/1.,  2.,  5.,  10./,  SIGN/1.,  -1./ 

OATA  NAHE/8H0*  T»£  3R98HA3  7,  UN,  8HTR EATED  ,8HC0TT0N  D, 

♦  dHUCK,  d.2»8H'.  (,l .  ,  3HTAN  46,  ,3HC0TT0N  P,8H0PLIN,  4, 

♦  8H  OZ.  ,8H  ,8H  , 8HT A  N  M-1,,8H  CL.  POL, 

♦6HY/WG0L  T , 8HR0P.  9  0,8HZ.  ♦  8H  ,8 HARM Y  GRE, 

♦8HEN  44,  W, 8HL  GAB.,  ,8H130Z.  , 8H  ,8H 

♦8H0LIVE  GR, SHEEN  107,, 8H  NYLON/C, 8HTN  (50/5, 8H0)  POPLI* 

♦  8HN,  5  OZ .  « 8H8LUE  150.8HTROP,  WL,  8H.  10  0Z.,8H  , 

♦8H  ,8H  , 8H ARMY  GRE, 8HEN  344,  ,8HP0LYESTE, 

♦  8HR/WOOL  G.8HA3AR0INE, dH  9  OZ.  , 3HARMY  GRE, 8 HEN  344,  , 

♦  8HPOL YESTE.8HR/ WOOL  T, 8HROP  ICAL, , 8H  10  OZ .  , 8HARM Y  GRE, 

♦8HEN  44,  W,8HL.  SERGE, 8H  15  OZ.  , 8H  ,8H  , 

♦8HOLIVE  GR.8HEEN  13d , , 3H  WO  CL/NY, 8HLON,  16  ,8H0Z.  SHIR, 

♦  3HTING  .dHCLIVE  GR,dHEEN  107, ,8H  CTN.  SA  ,  8HTEEN,  8., 

♦8H2  OZ.  ,8H  , 8H3LUE  151, 8H,  TROP.  ,8HWL.,  10  , 

♦  8HOZ .  ,8H  ,8H  / 

REAO  LAB  VALUES  OF  STANDARDS  (BUT  IN  THE  ORDER  ABL) 
REWIND  3 

DO  600  NS TO  =  1,  12 

IF  ( NSTD  .  NE  .  2)  NP  =  9 

IF  (NSTD  .  EQ  .  2)  NP  =  8 

REAO  (3,  1001) (CIE(3 ,  j),  CIE(1,  J) ,  CIE(2,  J),  J  =  1,  NP) 

ASSIGN  STARTING  VALUES  TO  THE  SEVEN  PARAMETERS 

HOMAX  =  0. 

VOMAX  =  0. 

NPM1  =  NP  -  1 
DO  150  J  *  1,  NPM1 
JP1  *  J  ♦  1 
DC  150  K  =  JP1,  NP 

HO  =  SORT  ( (CIE(  1 ,  J)  -  Cl  E  ( 1,  K))  **  2  ♦  CIE(2,  J) 

♦  -  C I£(  2,  K)  )  **  2) 

IF  (HO  .  LE  .  HDMAX)  GO  TO  130 
HOMAX  =  HO 
UHMAX  =  J 
KHMAX  =  K 

130  VO  =  A3S ( C IE ( 3,  j)  -  CIE(3,  K)> 

IF  (V0  .  LE  .  VOMAX)  GO  TO  150 
VOMAX  =  VO 
jVMAX  =  J 
KVHAX  =  K 
150  CONTINUE 

PAR ( 1 )  =  HOMAX  /  2. 


001000 
001010 
001020 
001030 
001040 
001050 
001060 
901070 
0010  80 
001085 
001090 
001100 
001110 
001120 
001130 
001140 
001150 
001160 
001170 
001180 
001190 
001200 
001210 
001220 
001230 
00124Q 
00125Q 
001260 
001270 
001280 
001290 
001300 
001310 
001320 
001330 
001340 
001345 
001350 
001360 
001370 
001364 
001390 
001400 
031410 
001420 
001430 
001435 
001440 
001445 
001450 
001460 
001470 
0015 01 
001510 
001520 
0  015  30 
00154Q 
001550 
001560 
001570 
001580 
001590 
001600 
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C  HORIZONTAL  MAJOR  SEMIAXIS 

PAR ( 2 )  =  PAR ( 1)  /  2. 

C  HORIZONTAL  MINOR  SEMIAXIS 

PAR(3>  =  VOMAX  /  2. 

C  VERTICAL  SEMIAXIS 

PAR ( 4 )  =  0.5  *  (CI£(1,  jHMAX)  ♦  CI£(i,  KHHAX) I 
C  A  STAR  COORDINATE  OF  CENTER  OF  ELLIPSOID 

pAR ( 5 i  =  0.5  *  (C IE t 2,  JHMAX)  f  CIE(2,  KHMAX)) 

C  8  STAR  COORDINATE  OF  CENTER  OF  ELLIPSOID 

PAR (6)  =  0.5  *  (CIE(3,  JVMAX)  ♦  CIE(3,  KVMAX)  ) 

C  L  STAR  COORDINATE  OF  CENTER  OF  ELLIPSOIO 

PAR ( 7 1  =  ATAN ((CIE(2 ,  JHMAX)  -  CIE{2,  KHHAX))  / 

♦  (CIE  (1  *  JHMAX)  -  Cl  £  ( 1 ,  KHMAX))) 

ANGLE  OF  HORIZONTAL  MAJOR  AXIS  TO  A  STAR  AXIS 

00  LOOP  FOR  IOIOT  SEARCH  PROCEDURE.  CHANGE  ONE 
PARAMETER  AT  A  TIME  ANO  SEE  IF  THE  ELLIPSOIO 
GETS  SMALLER.  IF  SO,  START  OVER.  EVENTUALLY 
MAKE  THE  PARAMETER  CHANGE  INCREMENTS  SMALLER 
AND  SMALLER. 

CALL  SIZER  (FAR,  ELL MIN) 

DO  450  NO  I V  -  1,  4 
L60  CHANGE  =  .FALSE. 

00  400  NPAR  =  1,  7 
DEL  =  DEL T A (NPAR)  /  OIV(NOIV) 

00  400  NSIGN  =  1,  2 
IF  (NSIGN  .  EQ  .  1)  GO  TO  170 
IF  (ACTION)  GO  TO  <,0  0 
170  ACTION  =  .FALSE. 

180  PAR  (NPAR)  =  PAR  (NPAR )  »•  OEL  *  SIGNCNSIGN) 

IF  ( PAR ( 1 )  .  LT  .  0.  .  OR  .  PAR(2)  .  LT  .  0 . 

♦  .  OR  .  PAR ( 3 )  .  LT  .  0.)  GO  TO  390 
CALL  SIZER  (FAR.  ELL Ic) 

IF  ( ELL  IE  .  GT  .  ELLMIN)  GO  TO  390 
ELLMIN  =  ELL  IE 
CHANGE  =  .TRUE. 

ACTION  =  .TRUE. 

GO  TO  180 

390  PA R ( NPAR)  =  PAR(NPAR)  -  OEL  *  SIGN(NSIGN) 

<.00  CONTINUE 

IF  (CHANGE)  GC  TO  160 
450  CONTINUE 

NORMALIZE  ELLIPSOIO  SEMIAXES  BEFORE  PRINTING  OUT 
INFORMATION.  SEMIAXES  ARE  AOJUSTEO  SO  AS  TO  MAKE 
RIGHT  HAND  SI OE  OF  ELLIPSOIO  EQUATION  UNITY. 

CALL  SIZER  (PAR.  ELL  IE ) 

ETEMP  =  PAR(l)  *  PAR ( 2 )  *  PAR(3) 

SCFACTR  =  (ELLIE  /  ETEMP)  **  (1.  /  3.) 

00  500  J  a  1.  3 
500  PAR ( J)  =  PAR ( j)  *  SCFACTR 

IF  (PAR(l)  .  GT  .  PAR{ 2) )  GO  TO  510 
TEMP  =  PAR (2) 

PA  R ( 2 )  =  PAR(l) 

PAR(l)  =  TEMP 

PAR ( 7 )  =  PAR ( 7)  -  0.5  *  3.141592654 

CALCULATE  ELLIPSOIO  PARAMETERS  FOR  OUTPUT 

510  ANG  =  P AR ( 7)  *  180.  /  3.  141592654 
IF  (ANG  .  GT  .  90)  ANG  =  ANG  -  180. 

SINE  =  SIN(PAR(7) ) 

COSINE  =  COS (PAR ( 7 ) ) 

Gil  =  (COSINE  /  PAR( 1) )  ••  2  ♦  (SINE  /  PAR(2)>  **  2 
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001610 
001620 
001630 
001640 
001650 
001660 
OQ 1670 
001680 
001690 
001700 
001710 
001720 
0017  30 
001740 
001750 
001760 
001770 
001780 
001790 
001800 
0018  30 
001840 
0018  60 
001865 
001870 
001880 
001890 
001900 
001910 
001920 
001930 
001935 
0019  36 
001940 
001950 
QO 1960 
001990 
00  20  00 
002010 
0020  20 
0020  30 
002040 
002050 
0020  60 
002070 
002080 
002090 
0  0  2100 
002110 
002120 
302130 
002140 
002150 
002153 
002156 
002159 
0021 62 
002165 
002168 
002170 
002180 
002190 
002195 
302200 
062210 
002220 


o  o  o  no  ooo  noo 


T HOG  12  =  2.  *  SINE  *  COSINE  *  Cl.  /  PAR(l)  **  2  -  002230 

♦  1.  /  PAR (2)  **  2)  002240 

G2  2  =  (SINE  /  PAR (1) )  *•  2  ♦  (COSINE  /  PARC  2) )  **  2  002258 

G33  =  l.  /  PARC  3)  **  2  002260 

VOL  =  <4.  /  3.)  *  3.  14  1592654  *  PARCll  ♦  PARC2)  *  PARC3)  002270 

002310 

PRINT  OUT  RESULTS  002320 


002330 

002340 

002350 

002360 

002370 

002380 

002385 

002386 

002307 

002388 

DEL  2  *  002390 

002395 

H0WFARC  J)  «  002400 

002403 

002405 

002406 

002406 


FORMAT  STATEMENTS  002409 

002410 

1001  FORMAT <  3F  8 . 2 )  002420 

1002  FORMAT (//////•  ACCEPT A3ILIT  T  ELLI  FSOI  0*  /IX,  6A8)  002430 

1003  FORMAT ( //24H  CENTER  OF  ELLIPSOIOl  L*,  F18.2/22X,  002440 

*■  2HA*,  F18.2/22X,  2H fl* »  F18.2)  002450 

1004  FORMATC/*  COEFFICIENTS;*, 8X,  *G11*.  E23.4/21X,  002460 

♦  *2G1 2*.  £2 3 • 4/ 22 X,  *G22*,  E23.4/22X.  *G33*.  E23.4I  002470 

1005  F ORM AT ( /*  SEMIMAJOR  HORIZONTAL  AXIS,  LENGTH*.  F8.2/  002480 

♦  *  SEMIMAjOR  HORIZONTAL  AXIS,  ANGLE*,  F8. 1/*  SE MI  MINOR* ,  002490 

♦  *  HORIZONTAL  AXIS,  LENGTH*,  F8.2/*  SEMIvERTICAL  AXIS,*  002500 

♦  *  LENGTH*,  F16.2)  002510 

10C6  FORMAT ( /♦  VOLUME  OF  ELLIPSO  I0*,F22. 2>  002520 

L  307  FGRMATC/*  OISTANCE  OF  POINT  FROM  CENTER  OF  ELLIPSOID*  0C2523 

♦  *  IN  ELL  UNITS*/ /IQ X,  *COOROINATES  OF  POINT  OISTANCE*,  002525 

♦  //11X,2HA*,6X,2H8*,6X,2HL*/C6X,  3F8.2,  Fli.5>)  002527 

LOO  8  FORMAT (////////)  002528 

end  0025  30 

002540 
0025  50 

SUBROUTINE  SIZER  (PAR,  ELL  IE)  002560 

002570 

DETERMINES  THE  SIZE  OF  THE  ELLIPSOID  PASSING  THROUGH  002580 

EACH  OF  THE  STANOAROS  WITH  THE  CURRENT  PARAMETERS,  AND  002590 

C  SELECTS  THE  ELLIPSOID  OF  LARGEST  SIZE  AMONG  THEM.  002600 

C  002610 

DIMENSION  PAR  <7 >  002620 

C  0  026  30 

COMMON  C I E ( 3 ,  9),  NP  002640 

C  002650 

SINE  s  SIN(PAR{7>  >  002660 

COSINE  =  COS (FAR ( 7) )  002678 

Gil  =  (COSINE  /  PARC1I)  •*  2  ♦  (SINE  /  PARC  2)  )  **  c  002680 

T HOG  1 2  =  2,  *  SINE  *  COSINE  *  Cl.  /  PARC1)  **  2  -  002690 

♦  1.  /  PARC2)  **  Z)  002700 

G22  =  (SINE  /  PAR  (1)  )  **  2  ♦  (COSINE  /  PAR<2>>  **  2  002710 

G33  =  1,  /  PAR(3>  ♦*  2  002720 

RHSMAX  =  0.  302730 

DO  150  J  =  1,  NP  002740 

DELTA  =  CIE(1,  J)  -  PA R ( 4 )  002750 

CELT  9  =  C  IE ( 2 ,  J)  -  PAR(5)  002760 

OELTL  =  CIE(3,  J)  -  PAR(6>  ,,  002778 


PRINT 

100 

2, 

(NAME 

(J,  NSTO)  , 

J 

=  1, 

6) 

PRINT 

100 

2, 

PAR  (6 

),  PARC 

4)  , 

PARCS) 

PRINT 

100 

4, 

Gil, 

THOG12, 

G  2 

2, 

G33 

PRINT 

100 

5, 

P  AR  ( 1 

),  ANG, 

PA 

R  ( 

2).  PAR ( 3) 

PRINT 

100 

6, 

VOL 

00  5  30 

J 

=  1 

,  NP 

OE  LI  = 

Cl 

E  (1 

.  J) 

-  P  AR  ( 4 

) 

JEL2  = 

Cl 

E  (2 

,  J) 

-  PAR  (5 

) 

DEL  3  = 

Cl 

E  (3 

,  J ) 

-  PAR (6 

) 

HO  V.  FAR 

( J) 

= 

G 11  * 

DELI  * 

*  2 

♦ 

TW0G12  * 

DELI  * 

G22  * 

OE 

L2 

**  2 

♦  G33  * 

OE 

L3 

**  2 

PRINT 

1  GO 

7, 

(Ciec 

1,  J)  , 

CIE 

C  2 

,  J )  , 

CIEC 

3,  J ) , 

♦  J  =  1,  NP) 
600  CONTINUE 
PRINT  1006 


ouc  -  Ml  *  n£LTA  **  2  ♦  TW0G12  *  OELTA 
.*&'•  KLTB •• t  ♦  *  OtLTt  2 

IF  ( RHS  •  LT  •  RHSNAX)  GO  TO  150 
RHSMAX  *  RHS 

150  1u!eN=ERHSH»X  ••  1.5  •  oariu  •  paru> 
return 
ENO 


*  DELT3  ♦ 


•  PAR<3» 


0  0  27  30 
0  0  27  90 
002600 
0026 10 
0026  20 
0026  30 
002640 
0026  50 


5- 
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PROGRAM  ELPS002  ( INPUT »  OUTPUT,  TAPE3) 

FOR  A  SERIES  OF  POINTS  IN  CIELA8  COLOR  SPACE, 
DETERMINES  THE  SMALLEST  NON-TILTING  ELLIPSOID 
THAT  CONTAINS  THE  POINTS.  THE  ELLIPSOID  MAY 
NOT  TILT  EITHER  FORWARD  OR  SIOEWAYS,  AND  MUST 
HAVE  ITS  MAJOR  AXIS  INCLINEQ  IN  THE  DIRECTION 
OF  THE  HUE  ANGLE  OF  THE  STANDARD. 

DIMENSION  PAR  (7  I  ,  HOWFARI9),  0IV(4),  NAME  <6,  121,  SIGN<2) 

COMMON  CIE(3,  91 ,  NP 

LOGICAL  CHANGE,  ACTION 

OATA  CELT  A/0* 1/,  OIV/1.,  2.,  5.,  10./,  SIGN/1.,  -1./ 

DATA  NAME/3HOLIVE  OR, 3H  A3  7,  UN.dHTREATEO  ,8HC0TT0N  D, 

♦  8HUCK,  8. 2 « 8H5  OZ .  ,  3HT  AN  46,  ,8HC0TT0N  P, 8HOPLIN,  4, 

*8H  OZ.  ,8N  ,  9H  , 8HTAN  M-1,,8H  CL.  POL* 

♦eHY/WCUL  T.8HROP.  9  0, 8HZ.  , 6H  »8HARMY  GRE, 

♦  3HEN  44,  W,  3HL  GAB.,  ,3H1T,0Z.  ,  3H  ,3H  , 

♦  8HOL I VE  GR,8HEEN  107,,  3H  NYLON/C, 8HTN  (50/5, 8H0)  POPLI, 

♦  8HN,  5  OZ . ,  8HELUE  15Q,8HTR0P.  HL,«H.  10  OZ..8H  , 

♦8H  ,8H  , 3HARMY  GRE.8HEN  344,  , iHPOLYESTE, 

fdHR/WOOL  G, JHA8ARQINE, 8H  9  OZ.  , 8HARMY  GRE, 8 HEN  344,  , 
♦8HPOLYESTE , 8HR/WOOL  T, 8HR0P ICAL, , 8H  10  OZ.  , SMARMY  GRE, 
♦8HEN  44,  W,3HL.  SERGE, 8H  15  OZ.  , 8H  ,8H  , 

8H0LIYE  GR,8HEEN  103,, 3H  WOOL/NY, 3HL0N,  16  ,3HOZ.  SHIR, 
♦8HTING  ,8H0LI/E  GR, 8HEEN  107,, 8H  CTN .  SA,8HT EEN,  8 ., 
♦8H2  OZ.  ,8H  , 8HBLUE  151, 6H,  TROP.  ,8HHL* ♦  11  , 

♦8HOZ.  ,8H  ,8H  / 

READ  LAB  VALUES  OF  STANDARDS  (BUT  IN  THE  ORDER  ABU 

REWIND  3 

DO  600  NSTD  =  1,  12 
IF  (NSTD  .  NE  .  2)  NP  =  9 
IF  (NSTO  .  EQ  .  21  NP  =  8 

READ  (3,  100 1) (Cl E(3 ,  J)  ,  CIE(1,  J) ,  CIE(2,  J),  J  =  1,  NP) 

ASSIGN  STARTING  VALUES  TO  THE  SEVEN  PARAMETERS 

PAR(4)  =  0. 

PAR ( 5)  =  0. 

PAR(6)  =  0. 

DO  120  J  =  1,  NP 

PAR ( 4)  =  PAR ( 4)  ♦  Cl E( 1,  Jl  /  NP 
C  A  STAR  COORDINATE  OF  CENTER  OF  ELLIPSOID 

PAR ( 5 )  =  PAR (5)  ¥  CIE( 2,  Jl  /  NP 
C  B  STAR  COORDINATE  OF  CENTER  OF  ELLPISOIO 

PAR ( 6 )  =  PAR (61  ¥  CIE( 3,  Jl  /  NP 
C  L  STAR  COORDINATE  OF  CENTER  OF  ELLIPSOID 

120  CONTINUE 
DiMAX  =  0. 

02MAX  =  9. 

03MAX  =  0. 

DO  150  J  =  1,  NP 

O  =  SORT ( ( PAR (4 1  -  CIE  (1 ,  Jl)  **  2  ♦  (PAR(5) 

♦  -  C  IF.(  2*  Jll  **  21 

ALPHA  «  ATAN  ( (PAR  (5)  -  CIE(2,  J)  I  /  (PAR(4> 

♦  -  CIE(  1,  Jl  1 1 

THETA  =  ATAN (PAR(5)  /  PAR(<»I| 


0010  00 
001010 
001020 
001030 
001040 
0  010  50 
0  01060 
001070 
001030 
001090 
001100 
001110 
001120 
001130 
001140 
001150 
001160 
001171 
001181 
001190 
011200 
001210 
301220 
001230 
001240 
001250 
001260 
001270 
001231 
001290 
001300 
081310 
101321 
001330 
•01340 
081351 
001360 
001370 
001380 
C 01398 
001400 
001410 
001420 
001430 
001440 
001450 
001460 
001470 
001483 
0  81490 
001500 
001510 
0015  20 
0015  30 
001540 
0015  50 
001560 
001570 
001581 
001590 
011600 
091610 
0  016  20 
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PHI  =  ALPHA  -  THETA  0016  30 

01  =  A3S ( 0  *  COS  I  PHI ) t  001640 

IP  101  .  LE  .  Q1MAX)  GO  TO  130  001650 

01  MAX  =  01  001660 

130  02  =  ASS ( 0  SI N ( PHI ) I  001670 

IP  (02  .  LE  02MAX)  GO  TO  140  001680 

G2MAX  =  02  001690 

l'*0  03  =  A8S  (  PAR  (6)  *  Cl  EC  3*  j»»  001700 

IP  (03  .  LE  .  03MAXI  GO  TO  150  001710 

J3MAX  =  03  001720 

150  CONTINUE  001730 

PAR(1»  =  01MAX  001740 

C  HORIZONTAL  MAJOR  SEMIAXIS  001750 

PAR(2I  =  02MAX  001760 

C  HORIZONTAL  MINOR  SEMIAXIS  001770 

PAR (  31  =  D3MAX  001780 

C  VERTICAL  SEMIAXIS  001790 

PAR(7)=  THETA  00  IS  00 

ANGLE  CP  HORIZONTAL  MAJAR  AXIS  TO  A  STAR  AXIS  001610 

001820 

00  LOOP  FOR  IDIOT  SEARCH  PROCEDURE.  CHANGE  CNE  001830 

PARAMETER  AT  A  TIME  AND  SEE  IF  THE  ELLIPSOID  001840 

GETS  SMALLER.  IF  SO,  START  OVER.  E/ENTUALLY  001850 

make  the  parameter  change  increments  smaller  ooiogo 

ANO  SMALLER.  001870 

001880 

CALL  SIZER  (PAR,  EllMIN)  001890 

CO  450  NOIZ  =  1»  4  '  001900 

160  CHANGE  =  .FALSE.  001910 

00  4 1: 0  NPAR  =  1,  6  001920 

DEL  =  DELTA  /  OlVCNOlV)  001930 

CO  400  NSIGN  =  1,  2  001940 

IF  (NSIGN  .  EQ  .  1)  GO  TO  170  001950 

IF  ( ACT  ION )  GC  TO  -,00  001960 

170  ACTION  =  .FALSE.  001970 

180  PAR  <  NPARl  =  PAR  (NPAR )  DEL  *  SIGN(NSIGN)  001980 

IF  (PAR(1»  .  LT  .  0.  .  OR  .  PARC2)  .  LT  .  0 .  001990 

♦.  OR  .  PARC3)  .  LT  .  Q.»  GO  TC  390  002000 

IF  (NPAR  .  EQ  .  4  .  OR  .  NPAR  .  EQ  .  5>  002010 

♦  PAR(7»  =  ATAN(PAR(5)  /  PAR  (4)  )  002020 

CALL  SIZER  (PAR,  EtLIE)  002030 

IP  (ELLIE  .  GT  .  ElLMIN)  GO  TO  390  002040 

EcLMlN  =  ELLIE  002050 

CHANGE  =  .TRUE.  002060 

ACTION  =  .TRUE.  002070 

GO  TO  18G  002080 

390  PAR  <  NPAR)  =  PAR  (NPAR)  -  DEL  *  SIGN(NSION)  002090 

IF  (  .PAR  .  EQ  .  4  .  OR  .  NPAR  .  E Q  .  5 »  002100 

*  PAR  ( 7 )  =  AT  AN  (  PA  R  (5  )  /  PAR(4J>  002110 

400  CONTINUE  002120 

IF  (CHANGE  >  GO  TO  loO  002130 

450  CONTINUE  002140 

002150 

NORMALIZE  ELLIPSOID  SEMIAXES  BEFORE  PRINTING  OUT  002160 

INFORMATION.  SEMIAXES  ARE  ADJUSTED  SO  AS  TO  MAKE  002170 

RIGHT  HANG  SIDE  OF  ELLIPSOID  EQUATION  UNITY.  002180 

002190 

call  SIZER  (PAR,  ELLIEI  002200 

fTEHP  =  PAR(l)  *  P..R<2»  *  PARC3)  002210 

SCFflCTR  =  (ELLIE  (  £TEMP»  **  Cl.  /  3.)  002220 

DO  530  j  =  1,  3  002230 

500  PAR ( J )  =  PARCJ)  •  3CFACTR  002248 

c  002250 

C  CALCULATE  ELLIPSOID  PARAMETERS  FOR  OUTPUT  002260 

c  002270 

ANG  =  P AR ( 7 1  *  180.  /  3.  14 1592654  002280 

!>4 


ooooooo  noo  ooo 


IF  (ANG  .  GT  .  90)  ANG  =  ANG  -  18C. 

SINE  =  SIN(PAfiC7) ) 

COSINE  =  COS (FAR (7) ) 

Gtt  =  (COSINE  /  PAR(l))  ♦*  2  ♦  (SINE  /  PAR(2)>  ♦*  2 
TWOG12  =  2.  *  SINE  *  COSINE  *  Cl.  /  PAR(l)  **  2  - 

♦  1.  /  PAR ( 2)  **  2) 

G22  =  (SINE  /  PAR  (1 )  )  **  2  ♦  (COSINE  /  PAR(2))  **  2 
G33  =  1.  /  PAR( 3)  **  2 
RATI01  =  PAR  ( 3)  /  PA R ( 2 ) 

RATI02  =  PAR (1)  /  PA  R ( 2 ) 

VOL  =  (4.  /  3.)  *  3,  141592654  *  PAR(l)  *PAR(2>  *  PAR(3> 
PRINT  OUT  RESULTS 

PRINT  ICQ  2,  ( NAME  (j,  NSTD),  J  =  1,  6) 

PRINT  100  3,  PAR  (6  ) ,  PARU),  PAR  <5 1 
PRINT  1004,  Gil,  TWOG12,  G22,  G33 
PRINT  1005,  PAR(1I,  ANG,  PAR(2I,  PAR(3) 

PRINT  1006,  RATIOl,  RATI02 

PRINT  1007,  VCL 

OO  53C  J  =  1,  NP 

CEL  1  =  CIE  (1 ,  J)  -  P AR (4 ) 

0EL2  =  Cl E (2  ,  J)  -  P  AR ( 5 ) 

DEL3  =  CIE (3 ,  J)  -  P AR (6 ) 

530  HOWFAMJI  =  Gil  *  DELI  *♦  2  ♦  THOG12  *  OELl  *  0EL2  * 

♦  G22  ♦  OEL2  **  2  ♦  G33  *  0EL3  2 

PRINT  1006,  (CIE (1,  J),  C I E (2 ,  J) ,  CI£(3,  J),  HOWFAR(J), 

♦  J  =  1,  NP) 

600  CONTINUE 

PRINT  1009 

FCKMAT  STATEMENTS 

1001  FORMAT ( 3F  8 . 2 ) 

1002  FORMAK//////*  ACCEPT  A  31  LI  T  Y  ELLI PSOI D*/l X, 6 A8> 

1003  FORMAT ( 7/ 24H  CENTER  OF  ELLIPSOID*  L*  »  F18.2/22X, 

*■  2HA * ,  F18.2/22X,  2h8*,  FJ8.2) 

1004  FORMAT!/*  COEFF IC  IENTS  **  ,  13  X  ,  *Glt*»  E23..4/21X, 

♦  *2G12* ,  E23.4/22X,  *5 22*,  E23.4/22X,  •G33*,  E23.4) 

1005  FORMAT ( /*  SEMIMAJOR  HORIZONTAL  AXIS,  LENGTH*,  F8.2/ 

♦  *  SEMIMAJOR  HORIZONTAL  AXIS,  ANGLE*,  Ftt.  1/*  SEMIMINOR*, 

♦  *  HORIZONTAL  AXIS,  LENGTH*,  Ffi.2/*  SEMIVERTICAL  AXIS,* 

♦  •  LENGTH*,  F16.2) 

1006  FORMAT!/*  LIGHTNESS  1  CHROMA  *  HUE  RATIOS*, 

♦  F10.2,  *  **,  F5.2,  *  *  1.00*) 

1007  FORMAT!/*  VOLUME  OF  ELLPISOIO*,  F22.2) 

1008  FORMAT!/*  DISTANCE  OF  POINT  FROM  CENTER  OF  ELLIPSOID* 

♦  *  IN  ELL  UN ITS*//1Q X,  ’COORDINATES  OF  POINT  DISTANCE*, 

♦  //11X,  2HA’,6X,  2H3*,  6X,2HL*/(6X,  3F8.2,  Fll.5>5 

1009  FORMAT! ////////) 

END 


SUBROUTINE  SIZER  (PAR,  ELLIE) 

DETERMINES  THE  SIZE  OF  THE  ELLIPSOIO  PASSING  THROUGH 
EACH  OF  THE  STANOARLS  HITH  THE  CURRENT  PARAMETERS,  AND 
SELECTS  THE  ELLIPSOID  OF  LARGEST  ...  ZE  AMONG  THEM. 

DIMENSION  PAR (7 ) 

C 

COMMON  CIE (3,  9) ,  NP 
C 

SINE  =  SIN (PAfi( 7 )  ) 

COSINE  =  COS ( PAR( 7) ) 

Gil  =  (COSINE  /  PAR(l))  *•  2  ♦  (SINE  /  PAR(2))  **  2 


002290 
002300 
002310 
0  0  2  3'.J 
002330 
002340 
002350 
002360 
002370 
002380 
002390 
002400 
002410 
002420 
002430 
002440 
002450 
002460 
00247* 
002480 
002490 
002500 
00251* 
0025  20 
0  025  30 
0  0  254* 
002550 
002560 
0025  70 
002580 
002590 
0026  00 
002610 
0026  20 
0026  33 
002640 
002650 
0  0  2660 
302670 

t/u  CE80 

0  C  2690 
002700 
002710 
002720 
0027  30 
002740 
0  0  27  50 
0  0276* 
0027  7* 
002780 
002790 
0  0  28  00 
00281* 
0  0  28  20 
002830 
002840 
0C2650 
0  0  28  60 
0  0  28  70 
002880 
00  26  91 
00  29*0 
8*291* 
0029  20 
0  0  29  39 

00294* 


TW0012  =  2.  *  SINE  •  COSINE  *  (1.  /  PARtll  Z  - 

♦  1.  /  PAR  (2)  •*  2) 

G22  =  (SINE  /  PAR(ll)  **  2  ♦  (COSINE  /  PARC  2) >  ** 
G33  =  1.  /  PAR( 31  **  2 
RHSHAX  =  0. 

00  150  J  -  1,  NP 

OELTA  =  CIEli.  J)  -  PAR ( 4 1 

OELTB  =  C IE ( 2 1  J)  -  PAR(5> 

CELTL  =  C I£( 3 «  J)  -  PAR(6) 

RHS  =  Gil  *  OELTA  ♦*  2  *  THOC-12  *  OELTA  *  OELTB  ♦ 

♦  G22  *  OELTB  **  2  *  G33  *  OELTL  **  2 
IF  (RHS  .  LT  .  RHSMAXH  GO  TO  150 
RHSHAX  =  RHS 

ELLIE  =  RHSHAX  ♦*  1.5  ♦  PAR(l)  *  PAR(2»  *  PAR(3> 

RETURN 

ENO 


2 


OC  2950 

0Q296Q 

0  52970 

002960 

002990 

0030  00 

00  30 10 

00  30  20 

0030  30 

003040 

003050 

093060 

003U70 

003080 

00  30  90 

003100 

003110 


Appendix  C.  Computer  program  tor  rearranging  file  in  numerical 
order  (Program  (REARRAN) 

PROGRAM  REARRAN  IINPUF*  OUTPUT*  TAPEi*  TAPE31 
DIMENSION  I ( 2GQ  I  *  AI2Q3).  3C200) 

REAL  L(20Q) 

LOGICAL  OK 
REMIND  3 
J  =  1 

100  READ  (3*  10011  U*  12*  LtJl*  Atjl*  B(J) 

1001  FORMAT ( 15 .  IX*  13,  F11.0,  F15.0.  F10.Q1 
IF  <11  .  EQ  .  999)  GO  TO  120 
I<J)  =  1000  *  II  *12 
J  =  J  ♦  i 
GO  TO  100 

120  NCOLORS  =  J  -  1 

NMIN1  =  NCOLORS  -  1 
ITER  =  0 
130  OK  =  .TRUE. 

ITER  =  ITER  ♦  1 
PRINT  1 00 A ,  ITER 
1004  FORMAT ( 1HS • 14) 

DO  160  J  =  l,  NMIN1 

IF  <I(J)  .  LE  .  IU  ♦  1>>  GO  TO  160 

OK  =  .FALSE. 

ITEMP  =  I(J) 

IU)  =  I(J  ♦  1) 

I(J  *  1>  =  ITEMP 
TEMP  s  KJ) 

L  ( J)  =  L ( J  ♦  1) 

LIJ  ♦  1)  =  TEMP 
TEMP  =  A ( J) 

A  <  J)  =  A  ( J  ♦  1) 

A ( J  ♦  1)  =  TEMP 
tem;-  =  b  <  j) 

3 ( J)  =  3( J  ♦  1) 

B(J  ♦  1)  =  TEMP 
160  CONTINUE 

IF  (.NOT  .  OK)  GO  TO  130 
DO  ldO  J  =  1,  NCOLORS 
H  =  If  J)  /  1C00 
12  =  I<J)  -  1000  *  II 

WRITE  <  1,  1002)  II.  12.  L  ( J )  *  A ( J )  .  8<Jl 
1002  FORMAT  ( 15 .  IH-,  13.  F11.2.  2F10.2) 

160  CONTINUE 
STOP 
END 


001000 
0010 10 
001020 
001030 
001040 
001050 
001060 
001070 
001066 
001090 
001100 
001110 
001120 
001130 
001140 
001150 
00 1160 
001170 
001180 
001190 
001200 
001210 
081220 
001230 
001240 
001250 
001260 
0C127O 
001200 
001290 
001300 
001310 
001320 
001330 
001340 
001350 
001360 
001370 
001360 
001390 
001400 
001410 
001420 
001430 
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Appendix  D.  Computer  program  for  selecting  sample  pairs  for  observers* 
evaluation  (Program  SIXWAYS) 

Following  is  the  derivation  of  the  equations  used  in  Program  SIXWAYS  to 
determine  if  a  given  sample  lay  no  more  than  0.1  CIELAB  unit  from  one  of  the  six 
iso  lines  radiating  from  a  given  standard. 

1.  Distance  between  a  point  and  an  iso  (hue- lightness)  line 

)|(  jj(  j|( 

Let  L  ,  a  ,  b  be  the  L  ,  a  ,  b  coordinates  of  the  standard.  The  iso 
s  s  s 

(hue-lightness)  line  is  the  straight  line  between  the  standard  point  and  the 

*  *  * 

achromatic  axis  at  the  L  lightness  level.  Let  L  ,  a  ,  b  be  the  L  ,  a  ,  b 

s  P  P  P 

coordinates  of  the  point.  The  asterisks  will  be  omitted  for  simplicity. 

First ,  raise  or  lower  the  point  vertically  until  it  has  the  same  lightness 

value  as  the  standard.  The  distance  moved  in  color  space  is  L  -  L  . 

b  P 

Now  drop  a  perpendicular  from  the  new  position  of  the  point  to  the  iso 
(hue-lightness)  line.  The  equation  of  the  iso  (hue-lightness)  line  in  the  constant 
L  plane  is 


b  =  (bs/as)a‘ 


(D-l) 


The  equation  of  the  perpendicular  from  the  new  position  of  the  point  is 


(b-  bp)/(a-  ap)  =  -(as/bg). 


(D-2) 


We  solve  Equations  D-l  and  D-2  simultaneously  to  get  the  point  of  intersection, 


which  we  will  call  (a^,  b.).  We  find  that 

“i  3  “s  (“p  as  +  bp  V  '  <as  +  ^  ' 


(D-3) 


bi=bs  <apas  +  bpbs>  ,(as  +  bs) 


(D-4) 


Now  the  distance  between  the  point  and  the  iso  (hue-lightness)  line, 


which  we  will  call  D^,  is 


Di  =  |(Ls  -  y 2  +  %  -  v2 + bP  -  v  j  • 
_ 


(D-5) 


Substituting  Equations  D-3  and  D-4  into  Equation  D-5,  we  obtain 

Dl=  <Ls- V2  +  (flsbp-apbs)2/(as  +  bs>  • 

2.  Distance  between  a  point  and  an  iso  (chroma- lightness)  line 


(D-G) 


We  first  raise  or  lower  the  point  vertically  until  it  has  the  same  lightness 
value  as  the  standard,  as  we  did  under  1. 

We  then  drop  a  perpendicular  from  the  point  to  the  iso  (hue-chroma)  line. 
The  equation  of  the  iso  (hue-chroma)  line  is 


(b  -  bg)  /  (a  -  as)  =  -(ag/bs). 


(D-7) 


The  equation  of  the  perpendicular  from  the  new  position  of  the  point  is 


(b  -  b)  /  (a  -  a  )  =  b  la  . 
p  p  s  s 


(0-8) 


We  solve  Equations  D-7  and  D-8  simultaneously  to  get  the  point  of  intersection, 
which  we  will  call  (a.,  b.).  We  find  that 

l  \ 


a.  =  a  (a2  +  b2)  +  b  (a  b  -  a  b  )  /  (a2  +  b2) 

l  s  s  s  sps  sp  /vs  s 


b.  =  b  (a2  +  b2)  -  a  (a  b  -  a  b)  (a2  +  b2) 
i  s  s  s  sps  s  P  '  s  s 


(0-9) 


(D-10) 


Substituting  Equations  D-9  and  D-10  into  Equation  D-5  (with  D^  substituted 
for  D^)  gives 


D„  -  (L  -  L  )2  +  (a2  +  b2  -  a  a  -  b  b  )2/ (a2  +  b2) 
2  s  p  s  s  p  s  p  s  /  v  s  s 


s  sps  p 


s  s 


D-  11) 


3.  Distance  between  a  point  and  an  iso  (hue-chroma)  line 


The  required  distance  is 


Ds=  <'\-y  +(bs 


-b  >2] 
P 


(0-12) 
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ooooo  ooo  ooo  oooo  oooo  ooo  oooo  oooooo 


r 


PROGRAM  SIXWAYS  (INPUT,  OUTPUT,  T ^P£i ,  TAPE7 ) 

IN  A  GROUF  OF  COLOR  POINTS  IN  CIELA8  SPACE,  DETERMINES 
ALL  COLORS  THAT  ARE  ON  ISOtHUE-LI GHTNESS) , 

ISO (HUE-CHROMA)  ANO  ISO (CHROMA-LIGHTNESS) 

LINES  RADIATING  FROM  EACH  POINT  IN  TURN. 

DIMENSION  a(400),  3(400),  NUMBERUOQ),  NUMBv40,  3,  2), 

♦  DE  (40,  3,  2),  NP(  3,  2) 

REAL  L( 40 0 ) 

READ  FILE  OF  COLOR  POINTS. 

J  =  1 

REWINO  1 
REWIND  7 

100  READ  (1,  1001)  NUMBER! J) ,  L(J),  A(J),  B(J> 

IF  ( NUMBER ( J )  .  EQ  .  7H9999999)  GO  TO  120 
J  =  J  ♦  1 
GO  TO  100 

120  NCOLOPS  =  J  -  1 

SET  UP  DO  LOOP  BASED  ON  EACH  POINT. 

CO  300  J  =  1,  NCOLORS 

PRINT  1014,  J 

WRITE  (7,  1002)  NUMBER (J ) 

SET  UP  OO  LOOP  BASED  ON  THE  THREE  ISO  LINES  IN 
COLOR  SPACE. 

DO  220  IV I A  *  1,  3 
HI  =  0 
M2  =  0 

SET  UP  OO  LOOP  BASED  ON  TESTING  EACH  POINT  TO  SEE 
IF  IT  IS  ON  THE  ISO  LINE. 

OO  200  K  =  l,  NCOLORS 
IF  (J  .  EQ  .  K)  GO  TO  200 
GO  TO  (130,  135,  140),  I/IA 

FORMULAS  FOR  DISTANCE  BETWEEN  COLOR  POINT  AND  ISU  LINE. 

130  OIST  =  SORT ( (L( J)  -  L( K) )  **24  (A(J)  *  B(K)  -  A  {  K)  *  B(J)) 

♦  **  2  /  ( A  ( J )  **2  4  8 ( J )  **  2)) 

GO  TO  150 

135  OIST  =  SORT  (  (L  ( J )  -  LOO)  **2  4  (A(J)  **  2  4  3(J>  **  2 

4  -  A  ( Kl  ♦  A  ( J)  -  B« )  *  B  ( J ) )  **  2  /  (A  ( J )  **2  4  3(J)  **  2)) 
GO  TO  150 

140  OIST  s  SORT  (  ( At  J)  -  A(O)  **2  4  <B(J>  -  3(K)>  **  2) 

IF  DISTANCE  IS  TOO  GREAT,  ELIMINATE  THE  POINT. 

150  IF  (OIST  .  GT  .  0.1)  GO  TO  200 

DETERMINE  ON  WHICH  SIDE  OF  THE  ISO  LINE  THE  POINT  LIES. 
CALCULATE  THE  COLOR  DIFFERENCE  BETWEEN  THE  POINT 
ANO  THE  REFERENCE  COLOR. 

GO  TO  (160,  165,  170),  IVIA 

160  B I  =  B(J)  *  (A(K)  ♦  A(  J)  4  B  (Kl  *  B(J))  /  (AO)  **  2  4  6(j) 


001000 
001010 
0010  20 
0010  30 
001040 
0  010  50 
001060 
0  010  70 
001080 
0 CIO  90 
001100 
001110 
301120 
001130 
001140 
001150 
001157 
001160 
001170 
0  0 11  BO 
001190 
001200 
001220 
001230 
001240 
001250 
001260 
001270 
001280 
001290 
001300 
001310 
001320 
001330 
001340 
001350 
001360 
001370 
001380 
001390 
001400 
0  01410 
001420 
001430 
001440 
001450 
001460 
001470 
001480 
001490 
001500 
0  015 10 
001520 
0  0 15  30 
001540 
001550 
001560 
001570 
0  01580 
001590 
001600 
001610 
001620 
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r>  o  O 


^  **  2) 

IF  (31  .  GT  .  3<J>)  GO  TO  180 
GO  TO  190 

165  AX  =  AC  J>  ♦  8  CJ I  *  <  A  C  K  >  *  BUI  -  A(.J>  *  B(KII  / 

♦  *  BUI  **  2 1 

IF  <  A I  .  GT  .  A  ( J  )  I  GO  TO  190 
GO  TO  180 

170  IF  (UK)  .  GT  •  L(J)I  GO  TO  180 
GO  TO  190 
180  Ml  *  fl  M 

NUMB (Ml t  I VI  A*  1)  =  NUMBER ( K) 

OE  (Ml*  UIA,  II  *  5QRT  ( (L  ( Jl  -  L(KI1  **  2  ♦  (A(Jt 

♦  **  2  ♦  (B(JI  -  B  Ml  I  **  21 
GO  TO  200 

190  M2  =  M2  «■  .1 

NUMB ( M2 1  I  /I A  t  2)  =  NUMBER ( Kl 

OE  ( N2 »  I V  I A  t  2)  =  SORT  ( ( L  (  J I  -  L(Kll  **  2  *  (A(JI 

♦  2  ♦  (B(Jl  -  a«l  I  **  21 


(AUI  ** 


-  A (<  1  I 


-  A  <  K 1 1 


200  CON. INUE 

NPd/IAt  II  -  Ml 
NP ( I  V I A  t  21  =  H2 
220  CONTINUE 


2 


PRINT  OUT  RESULTS  FOR  EACH  STANQARQ  POINT. 


NM AX  =  MAX0 ( NP( It  lit  NPIlt  21 t  NP(2.  lit  NP(2t  2lt  NP(3, 
♦  NP(3,  2) ) 


IF  (NMAX  .  EQ  . 
WRITE  (7,  10111 
OO  250  K  =  It  NMAX 
WRITE  (7,  10041 
CO  250  Ml  =  1,  3 
00  250  M2  s  It  2 
IF  (K  .  GT  .  NP (M 1 
IF  (Ml  .  EQ  .  I 
«•  NUMB (K »  Hit  M2) 


0)  GO  TO  270 


IF  (Ml  .  EQ  .  1 

♦  NUMB (K t  Ml,  M2) 
IF  (Mi  .  EQ  .  2 

♦  NUMB (K,  Ml,  M2) 
IF  (Ml  .  EQ  .  2 

♦  NUMB  (K  t  Ml,  M2) 
IF  (Ml  .  EQ  .  3 

NUM8  (K  t  Ml  »M2  ) 
IF  (Ml  .  EQ  *  3 

♦  NUM 3 (K ?  Ml  ,M2 1 
250  CONTINUE 

GO  TO  300 

270  WRITE  (7,  10121 
300  CONTINUE 
STOP 


t  M2)  I 
ANO  . 
OE  (K, 

and  . 

OE  (<  , 
ANO  . 
OE  (K , 
ANO  . 
Q£  <<  ♦ 
ANO  . 
OE  ( K  , 
ANO  . 
OE  (K, 


GO  TO  250 
M2  .  EQ  . 
Ml,  M2) 

M2  .  EQ  . 
Ml,  M2I 
M2  •  EQ  . 
Ml,  M2I 
M2  .  EQ  . 
Ml  ,M2) 

M2  .  EQ  . 
Ml  ,M2 ) 

MZ  •  EQ  . 
Ml,  M2) 


II  WRITE  (7,  1005) 
21  WRITE  (7,  1006) 
1)  WRITE  (7,  1007) 
21  WRITE  (7,  10081 
II  WRITE  (7,  10091 
21  WRITE  (7,  10101 


C 

c  FORMAT  STATEMENTS 

C 

1001  F0RMAT(2X,A7«F12.2,2F13.2) 

1002  FORMAT (////*  SAMPLE  NUMBER  *,  A7I 

1003  FORMAT (15  I 
1009  FORMATS  *) 

1005  FORMAT ( 1H ♦  ,  A7.F6.2) 

1006  F0RMAT(1H»,17X,A7,F6.2) 


1007  FORMAT ( 1H*» 34X, A7 ,F6.2I 

1008  F0RMAT(lH»,51X,A7tF6.2l 

1009  F0RMAT(lH*,6dX,A7,F6.2l 

1010  FORMAT(lH*,85X,A7,F6.2l 

1011  FORMAT ( *  PLUS  CHROMA* ,  5X,  *MINUS  CHROMA*,  7X, 
♦*PLU$  HUE*,  9X,  *MINUS  HUE*,  7X,  *PIUS  LIGHT*,  7X, 
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1), 


0016  30 
00164C 
001658 
001660 
081670 
001688 
001690 
001708 
0017 10 
001720 
0017  30 
001740 
001750 
001760 
001770 
0017  80 
001790 
001830 
0018  10 
001820 
0018  30 
001840 
001B50 
001860 
001870 
001888 
001890 
001900 
001905 
001910 
001920 
001930 
001940 
001950 
001960 
001970 
001980 
001990 
0020  00 
0020  10 
0020  20 
0020  30 
00  20  40 
0020  50 
002060 
002070 
002080 
002090 
002100 
002110 
002120 
002130 
002140 
002158 
002160 
002170 
002200 
002210 
002220 
002240 
002250 
002260 
002270 
002280 
002290 
002300 


♦•MINUS  LIGHT*)  002305 

1012  FORMAT C *  NO  POINTS  FOUND  IN  ANT  OIRECTION*)  002310 

1014  FORMATC 1H < »1 5 )  002315 

C  002320 

ENO  002330 
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Typical  portion  of  output  from  Program  Sideways 


SAMPLE  NUMBER  519-152 
PLUS  CHROMA  MINUS  CHROMA 

513-  10  .21 

524-610  .09 

967-160  .06 


PLUS  HUE  MINUS  HUE 

524-610  .09  967-160  .06 

5  28-662  .16 

536-835  .12 


PLUS  LIGHT  MINOS  LIGHT 

519-336  .66'  528-610  .09 

526-753  2.08  967-160  .16 

528- 576  .15 

529- 757  .29 

529-765  .31 

529-779  .61 


SAMPLE  NUMBER  519-176 

PLUS  CHRCMA  MINUS  CHROMA  PLUS  HUE 

513-  3C  .13  526-  60  .66  513-  30  .13 

26-135  .11  526-123  .11  526-135  .11 

28-610  .18  524-546  .10  528-586  .10 

531-666  .29 


MINUS  HUE  PLUS  LIGHT  MINUS  LIGHT 

526-108  .11  51-989*00.20  526-123  .11 

511- 167  1.01  52*-135  .11 

512- 696  .76 

513-  9  .19 

513-  10  .17 

526-  71  .56 

526-764  1.66 

526-768  1.68 

526-753  1.69 

528- 607  .75 

524-636  .42 

529- 776  .76 

530- 365  .47 


SAMPLE  NUMBER  519-223 

FLUS  CHROMA  MINUS  CHROMA  PLUS  HUE 

513-  5  .23  529-764  .05 

513-  10  .29 

526-165  .66 

524-116  .20 

524-576  .10 

529- 768  .05 

530- 351  .61 

536-980  .68 

536-436  .26 

536-835  .10 

536-929  .67 


MINUS  HUE  PLUS  LIGHT  MINUS  LIGHT 

528-576  .10  526-753  1.96  536-835  .10 

536-635  .10  528-576  .10 

529-768  .05 

529-757  .16 

529-765  .19 

529-779  .69 


o  o  o 


Anoendix  E.  Computer  program  for  randomizing  order  of  presentation 
cf  sample  pairs  (Program  FUNFER) 

PROGRAM  FUNFER  (INPUT*  OUTPUT,  TAPED 

PREPARES  RANUIMIZEO  LIST  OF  COMPARISONS  FOR  THE  INSPECTOR 

DIMENSION  NAME  (2,  24),  ITAL(24) 

LOGICAL  DONE 

OATA  NAME/ 

♦1001, 84 44, 84 44, 1001*7030, 74 3C ,7430,70 30, 

♦  7291, 84 43, 8443, 7251,7052,  700G , 70  CO ,7052, 

♦  0709, 70 23,7023,0709,0  947 ,3979, 0979, 08 47, 

♦0 €6 6, 06 50, 0650, 0866, 0006, 0837, 08 37 ,0006, 

♦0681,0984,0984,0681, 0714,0644, 0644,0714, 
♦5000,5327.5327,5000,1000,5310,5310,1000/ 

DO  ICO  1  =  1,  24 
100  ITALU*  -  0 

150  K  =  IFI  X  ( RANF  ( 0,  0 )  *  24.  ♦  1.) 

IF  (ITAL(K)  .  GE  .  t  a »  GO  TO  150 
WRITE  (i»  1001)  (NAME ( J,  X),  J  =  1,  2) 

ITAL(K)  -  ITAL(K)  ♦  1 
CONE  =  .TRUE. 

DO  200  J  =  1,  24 

IF  (ITALCJ)  .  LT  .  10)  DONE  *  .FALSE. 

200  CONTINUE 

IF  (OONE)  STOP 
GO  TO  150 

1001  FCRWAT(//1X,  14,  19) 

END 


001GGC 
001010 
3010  20 
J010  30 
301040 
001350 
0C136C 
001370 
301075 
001080 
001090 
001100 
•301110 
001123 
301130 
•301220 
3C1230 
001240 
301250 
301260 
301270 
001200 
031293 
001300 
001310 
001320 
001330 
001340 
301250 
00137G 


M4NNYLW  000081  LINES  PRINTED 
MANNYLW  300081  LINES  PRINTED 


///  END  OF  LIST  /// 
///  END  OF  LIST  /// 


LO  22 
LG  22 


******  17,25.57 
******  17.25.57 


Appendix  F.  Computer  program  for  calculation  of  50%  acceptability  limits 
by  logistic  function,  (Program  LOGIT). 


The  logistic  function  was  used  by  Berkson  (4,  5)  to  establish  the  dose 
of  a  drug  that  is  lethal  to  50%  of  the  population  exposed,  the  so-called  L.D.50 
value.  The  use  in  our  work  is  analogous:  we  wish  to  determine  the  color  difference 
at  which  a  given  inspector  will  pass  examples  50%  of  the  time.  The  papers  of 
Berkson  should  be  consulted  for  details  of  the  function  and  its  use;  we  will  just 
present  the  equations  here. 

The  logistic  function  is 


P  =  1  -  Q 


-(a  +  3  x) 


1  +  e 


(F-l) 


In  this  equation,  x  is  an  independent  variable.  In  Berkson's  work  with  drugs,  it 
is  the  logarithm  of  the  dosage;  in  our  work,  following  Berkson,  we  used  the 
logarithm  of  the  color  difference.  The  symbol  P  represents,  in  Berkson's  work, 
the  calculated  fraction  of  test  animals  killed  by  the  dose  corresponding  to  x;  in 
our  work  it  is  the  calculated  fraction  of  times  the  inspector  passes  the  pair  of 
samples  having  a  color  difference  corresponding  to  x.  The  quantities  a  and  3  are 
parameters . 

Equation  F-l  can  be  linearized  by  using  the  function  In  (P/Q),  which  is 
equal  to  a  +  3  *.  This  function  is  known  as  the  logit  of  P.  Thus,  plotting  the  logit 
of  P  against  x  should  yield  a  straight  line,  if  tne  assumptions  on  which  the  logistic 
function  is  based  are  valid. 


4  J.  Berkson,  "Application  of  the  Logistic.  Function  to  Bio-Assay,"  Am.  Stat. 

Assn.  J.,  39,  357-365  (1944). 

5.  J.  Berkson,  "A  Statistically  Precise  and  Relatively  Simple  Method  of  Estimating 
the  Bio-Assay  with  Quantal  Response,  Based  on  the  Logistic  Function,"  Am.  Stat. 
Assn.  J.,  48,  565-599  (1953). 


The  following  analysis  applies  to  one  inspector  looking  at  a  series  of  four 
color  differences,  where  each  color  difference  is  shown  to  him  10  times.  Let  p^ 
be  the  observed  fraction  of  passes  at  the  color  difference  corresponding  to  x^ 
(remember  that  P  was  the  theoretical  fraction  of  passes) ,  q.  =  1  -  p. ,  and 
1.  =  In  (p. /q^).  It  can  then  be  shown  that  good  estimates  of  a  and  8  (denoted  by 
a  and  b)  can  be  obtained  by  the  following  normal  equations,  where  the  summations 
are  taken  over  the  four  color  differences: 


o 

_  Epql  £pqx  -  £pqlx  £pqx  , 

a  ~  2  2 
£pq  £pqx'  -  (£pqx) 


(F-2) 


b  -  SB.ailpqfr  -  sgqxJEfli- 

Epq  £pqx2  -  (Epqx)2  (F‘3) 

Once  a  and  b  have  been  calculated ,  the  x  value  corresponding  to  a  50%  pass 
is  given  very  simply.  Since  In  (P/Q)  =  a+  Bx,ifP  =  Q  =  0.5  then  In  (P/Q)  = 

0  and  xcn  =  -(a/B).  We  thus  have 

DU 

A  E-n  =  Q"(a/b)  (F-4) 

DU 

since  x  represented  the  logarithm  of  the  color  difference.  The  symbol  AE^  refers 

to  the  color  difference  for  which  the  inspector  will  pass  the  sample  50%  of  the  time. 

The  standard  deviation  of  the  AE_n  value  is  calculated  by  the  following 

ou 

equations : 


06 


s l  ,  =  1/(10  l  pq), 

H 


(F-5) 


s 


2 

b 


1/ 


E 


10  Z  pq  (x  -  x) 


where  x  is  the  x  value  averaged  over  the  four  color  differences. 


(F-6) 


(F-7) 


sAE50  =  AE50sx50 


(F-8) 


where  sAr,  is  the  desired  standard  deviation. 

AE50 

In  entering  the  observers*  data  for  fraction  of  passes,  we  avoided  entering 
1  for  100%  passes  or  0  for  0%  passes ,  since  the  corresponding  logits  would  have  been 
infinite.  Following  Berkson,  we  used  0.95  and  0.05,  respectively. 
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PROGRAM  LOGIT  (INPUT*  OUTPUT,  TAPED 

COMPUTES  50  PER  CENT  ELLIPS0I0  BOUNDARIES 
BY  THE  METHOD  OF  LOGISTIC  FUNCTIONS. 

OIMENSICN  0E  (A*  6,  21,  Q(4),  IQ<4,  3),  W(4>,  050(3), 
♦ST0EVI3I 

REAL  LOGIT ( 4 ) ,  L0G0E(4) 

INTEGER  CCHHLL,  GUY,  SIZE,  COLOR 

DATA  0E/.4,  .6,  1*13,  1.52,  *4,  .0,  1.13,  1.52,  .4,  .77, 
♦1.17,  1.32,  .4,  .77,  1.17,  1.45,  .6,  1.17,  1.02,  2.29, 

♦.6,  1.17.  1.82,  2.29,  .5,  .92,  1.33,  1.92,  .5,  .92,  1.33, 
♦1.92,  .25,  .47,  .69,  1.01,  .25^  .40,  .69,  1.01,  .66,  1.32, 
♦1.95,  2.65,  .66,  1.32,  1.95,  2.65/ 

00  LOOPS  FOR  COLOR  STANDARO,  OIRECTION  IN 
COLOR  SPACE,  ANO  INSPECTOR. 

REMINO  1 

00  300  COLOR  *  1,  2 
IF  (COLOR  .  EQ  .  1)  PRINT  1003 
IF  (COLOR  .  EQ  .  2)  PRINT  1004 
PRINT  1012 

00  300  CCHHLL  »  1,  6 

GO  TO  (101,  102,  103,  104,  105,  106),  CCHHLL 

101  PRINT  1005 
GO  TO  110 

102  PRINT  1006 
GO  TO  110 

103  PRINT  1007 
GO  TO  110 

104  PRINT  1008 
GO  TO  110 

105  PRINT  1009 
GO  TO  110 

106  PRINT  1010 
110  AVX  *  0. 

00  120  SIZE  *  1,  4 

LOGDE  (SIZE )  *  ALOG(OE(  SIZE ,  CCHHLL,  COLOR)) 

120  AVX  =  AVX  ♦  LCGOE  (SIZE) 

AVX  *  0.25  *  AVX 
00  250  GUY  *  1,  3 

REAO  PERCENT  PASSES. 

READ  (1,  1002)  (Q  (SIZE) ,  SIZE  a  1,  4) 

00  130  SIZE  *  1,  4 

IQISIZE,  GUY)  a  100.  *  Q(SIZE)  ♦  .0001 
IF  ( IQ( SI ZE,  GUY)  .  GT  .  90)  IQISIZE,  GUY)  a  100 
IF  ( IOCSIZE,  GUY)  .  LT  .  10)  IQISIZE,  GUY)  a  0 
130  CONTINUE 

CALCULATION  OF  LOGISTIC  FUNCTIONS  ANO  HEIGHTS • 

00  150  SIZE  *  1,  4 
P*l.-  0  (SI ZE) 

M(SIZE)  *  P  *  QISIZE) 

150  LOGZT (SIZE)  a  ALOGIP  t  Q(SIZE)) 

C 

C  CALCULATION  OF  SUMMATION  VALUES. 


001000 
001010 
0010  20 
001030 
001040 
001050 
001060 
001070 
0  010  80 
001090 
001100 
001110 
001120 
001130 
001140 
001150 
001160 
001170 

ooiiao 

001190 
001200 
001205 
001210 
001220 
001230 
001240 
0  C 1250 
001260 
001270  , 
001280 
001290 
001300 
001310 
001320 
001330 
001340 
001350 
001360 
001370 
001380 
001390 
001400 
0  01410 
0014  20 
0014  30 
001440 
001450 
001460 
001470 
001488 
001490 
001500 
001505 
0  01510 
0015  20 
001530 
0  01540 
001550 
001560 
001570 
001580 
001590 
001600 


OOO  OOO  OOO  ooo  ooo 


c 


200 


SUM*  =  0. 
SUMKL  =  0. 
SUHWLX  =  0. 

sunwx  =  e. 

SUHHX2  a  Q • 


SUMWXX2  =  0. 

oo  200  size  =  i,  u 

SUHW  =  SUHW  ♦  W (SIZE ) 


SUHWL  = 
SUHWLX 
SUHWX  = 
SUMHX2 
SUMWXX2 


SUMWL  ♦ 
=  SUHWLX 
SUHW X  ♦ 
=  SUHWX2 


W(SIZE)  * 

*  W( SIZE) 
W (31 ZE )  * 

♦  W(SIZE) 


LC6IT (SIZE) 

*  LOGIT (SIZE) 
LOGOE(SIZE) 

*  L0G0E (SIZE) 


*  LOGOECSIZE) 


**  2 


=  SUMWXX2  ♦  W(SIZE)  *  (LOGOECSIZE)  -  AVX)  ♦♦  2 


CALCULATION  OF  SLOPE  AND  INTERCEPT  OF  LEAST  SOUARES  LINE. 

OENOM  =  SUMW  *  SUMHX2  -  SUHWX  **  2 
A  =  (SUHWL  *  SUHWX2  -  SUHWLX  *  SUHWX)  /  OENON 
0  =  (SUHW  *  SUHWLX  -  SUHWX  *  SUHWL)  /  OENOH 

50  PER  CENT  POINT. 

X50  =  -A  /  8 

05  0 ( GUY )  =  EXP(X5fl) 


STANOARO  CEVIATION  OF  59  PER  CENT  POINT. 


250 


VARAPR  =  1.  ✓  (10.  *  SUHW) 
l/ARB  at.  /  (10.  *  SUHWXX2) 
/ARX5Q  =  (t.  /  8  **  2)  *  UARAPR 
VAR050  =  VARX50  '  05  0 ( GUY)  **  2 
STOEVCGUYJ  =  SORT  ( VA R050 ) 
CONTINUE 


♦  VARB  •  ( X59  -  A/X)  ••  2) 


PRINT  RESULTS. 

PRINT  1011*  (OE (SIZE* 

♦GUY  =  1,  3),  SIZE  =  1* 
PRINT  1001,  OE ( 3 «  CCHHLL. 
♦  05  0 ( GUY ) ,  STOEV (GUY) »  GUY 
PRINT  1011,  QE(4*  CCHHLL* 
♦GUY  =  1,  3) 

300  CONTINUE 
STOP 


CCHHLL* 

2) 


COLOR),  (IQ (SIZE,  GUY) 


COLOR)*  (IQ  (3*  GUY) 
=  1,  3) 

COLOR) * (10(4*  GUY), 


1001 


1002 

1003 

1004 

1005 

1006 
1007 
1006 

1009 

1010 
1011 
10 1 2 


FCRHAT  STATEMENTS. 

FORMAT (F5 . 2*  16,  F7.2,  2H  ♦/lH*,  17X,  2H 

♦  15,  F7.2,  2H  ♦/  1H ♦ ,  36X,  2H  «♦,  F5.2,  15* 

♦  IH ♦  *  55 X,  2H  i*,  F5.2) 

FORMAT(F15.Q) 

FORMAT (////22X,  *  C  OL  OR  *  OLIVE  GREEN*) 
FORMAT ( ////26X,  *COLORl  TAN*) 

FORMAT ( //6X»  *PLUS  CHROMA*/) 

♦MINUS  CHROMA*/) 

♦PLUS  HUE*/) 

♦MINUS  HJE*/) 

♦PLUS  LIGHT*/) 

•MINUS  LIGHT*/) 

16,  2I19/F5.2,  16, 

DELTA*,  7X ,  *  8 1 LL  *  * 


-,  F5.2, 
F7.2,  2H 


♦  / 


FORMAT ( //6X, 

FORMAT ( //6X» 

FORMAT! // 6X, 

FORMAT (//6X, 

FORM AT ( // 6  X » 

FORMAT ( F5 , 2* 

FORMAT!///* 

♦  /3X,  *E''/4X*  3C4X,  *PCT.*,  3X, 
♦2X,  *a  ST0.  G£V •* ) ) 


2119) 

16X,  *MOE*,  15X,  *CAROL* 
♦BOUNDARY* )/6X,  3(2X,  *PASS*, 


END 
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001610 
0  016  20 
001630 
QQ1640 
001650 
001660 
001670 
001660 
001690 
0017P0 
0  017  10 
001720 
001730 
001740 
001750 
001760 
001770 
0  0 17  60 
001790 
001800 
0  01610 
001620 
001630 
001640 
0016  50 
001660 
001670 
001660 
001690 
001906 
001910 
001920 
001930 
0  01940 
001950 
001960 
001970 
001960 
001990 
002000 
002005 
002010 
002020 
0020  30 
002040 
002050 
002060 
002070 
0020  60 
002090 
002095 
002100 
002110 
0  02120 
0  02130 
0  0  2140 
002150 
002160 
002170 
002160 
002190 
002200 
002210 
002220 
002230 
902240 


Appendix  G.  Derivation  and  Meaning  of  the  Acceptability  Equation 

If  the  acceptability  equation,  Equation  1,  were  to  have  the  square  root 
bracket  removed,  and  if  A  A  were  set  equal  to  unity,  it  would  become  an  ellipsoid 
equation  in  which  the  center  of  the  ellipsoid  is  placed  at  the  standard  point  in 
CIELAB  space,  the  three  semiaxes  are  made  equal  in  size  to  the  tolerance  values 
determined  by  the  inspectors ,  and  the  angle  6  is  made  equal  to  the  hue  angle  of 
the  standard.  Thus,  in  Equation  B-3,  c,  the  chroma  tolerance,  replaces  p;  h, 
the  hue  tolerance  replaces  q;  and  v  ,  the  lightness  tolerance,  replaces  r.  The 
coordinates  of  the  standard,  a^bg  and  Lg,  replace  those  of  the  center  of  the 

)|i  j(i  ]|( 

ellipsoid,  ac>  bc  and  Lq.  With  these  changes,  Equations  B-4  to  B-8  become 
Equations  2  to  6,  and  Equation  1  modified  as  just  indicated  is  equivalent  to  Equation 
B-9.  This  equation  defines  what  we  call  the  "acceptability  ellipsoid." 

Now  if  we  substitute  the  coordinates  of  any  point  into  this  equation ,  we 

generally  will  no  longer  calculate  unity  but  some  other  value  that  we  have  called 
2  2 

S  in  Appendix  B  and  we  will  call  (A A)  here.  We  showed  in  Appendix  B  that  we 
have  defined  a  new  ellipsoid  that  passes  through  the  sample  point  in  which  the 
linear  distances  are  expanded  or  contracted  by  a  factor  of  S  (or  AA).  We  thus  see 
the  rationale  of  taking  the  square  root  of  the  expression  in  Equation  1;  it  is  that  we 
wish  to  calculate  the  factor  (AA)  by  which  the  linear  distance  from  the  standard  to 
the  sample  in  question  is  greater  or  less  than  the  linear  distance  along  the  seme 
line  but  from  the  standard  to  the  acceptability  ellipsoid.  If  A  A  is  greater  than  1, 
the  point  lies  outside  of  the  acceptability  ellipsoid;  if  it  is  less  than  1,  the  point 
is  inside  the  acceptability  ellipsoid. 


